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INTRODUCTION 


Finite  difference  methods  are  today  one  of  the  important  tools  for 
approximating  the  solutions  of  time-dependent  problems  governed  by  systems  of 
partial  differential,  equations .  This  is  ture  for  example  for  a  wide  spectrum 
of  pure  initial  value  and  mixed  initial -boundary  value  problems  in  the  field 
of  fluid-dynamics.  The  widespread  use  of  finite  difference  schemes  to  solve  such 
problems  has  increased  rapidly  since  the  early  1950' s,  in  response  to  the 
increased  capabilities  of  the  electronic  computers  which  execute  the  vast  amount 
of  calculations  needed  in  the  applications. 

Since  that  time  there  have  also  been  extensive  developments  in  the  analysis 
of  finite  difference  schemes.  A  prime  example  of  this  has  been  the  maturation 
of  the  mathematical  theory  needed  to  handle  numerical  approximations  to  linear 
systems  of  initial  value  problems.  The  concept  of  stability  of  a  difference 
scheme,  expressing  a  continuous  dependence  of  the  scheme-solution  on  its  initial 
values,  plays  a  major  role  in  the  above  mentioned  theory.  The  centrality  of 
this  concept  follows  from  the  Lax  equivalence  theorem  (see  for  example  [13, 

Chapter  3]  which  assures  the  convergence  of  a  numerical  computation  carried  out 
by  a  stable  finite  difference  scheme  consistent  with  a  well-posed  initial  value 
problem. 

Besides  the  pure  initial  value  problems  mentioned  above fwe  are  most 
interested  in  approximating  the  solution  of  mixed  Initial -boundary  value  problems, 
where  the  numerical  approximation  must  include  a  special  boundary  treatment  to 
fulfi]  the  boundary  conditions  imposed  on  the  problem.  Furthermore,  since 


i 
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practically  all  the  various  schemes  are  solved  in  a  finite  domain  of  the  grid, 
it  follows  that  numerical  boundary  conditions  must  be  added  also  in  the  case  of 
pure  initial  value  problems  so  that  the  solution  can  be  uniquely  determined.  It 
follows  that  in  all  cases  (pure  initial  value  and  mixed  initial -boundary  value 
problems),  the  overall  approximation  is  composed  of  a  basic  scheme  applied  at 
inner  grid  points  and  a  (different)  additional  algorithm  which  is  applied  locally 
at  the  boundary.  The  additional  boundary  treatment  which  determines  uniquely 
the  scheme-solution  is  sometimes  an  artificial  one  and  does  not  necessarily 
reflect  the  boundary  conditions  (if  any)  of  the  original  differential  system, 
so  it  may  cause  an  instability.  Indeed  it  is  known  (see  for  example  [13,  Chapter 
6]  [11,  Chapter  17])  that  even  if  the  basic  scheme  is  stable,  a  careless 

| 

numerical  boundary  treatment  may  render  the  total  computation  unstable. 

These  considerations  lead  us  into  the  area  of  the  stability  analysis  of 
approximations  to  hyperbolic  initial-boundary  value  problems.  One  of  the  most 

important  contributions  in  that  area  which  will  serve  us  as  a  general  reference  | 

on  the  subject,  is  the  1972  paper  by  Gustafsson,  Kreiss  and  SundstrBm  [6]  j 

which  is  a  generalization  of  an  earlier  paper  by  Kreiss  (1968)  [8],  The  analysis  j 

in  the  1972  paper  rests  on  a  new  stability  definition  (Definitions  3.2  and  3.3  in  ^ 

[6])  which  like  the  stability  definition  for  approximations  to  pure  initial 
value  problems,  is  obtained  by  a  discretization  of  a  corresponding  well-posedness 
condition  of  the  original  differential  equation  (see  for  example  [93).  This  new 
stability  condition  reflects  the  influence  that  the  boundary  values  have  on  the 
numerical  solution,  and  as  in  the  cue  of  pure  initial  value  problems,  it  serves  as 
a  sufficient  condition  for  the  convergence  of  a  (compatible)  consistent 
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approximation  [5].  The  main  result  of  Gustafsson  et.  al.  [6]  provides  an 
algebraic  criterion  which  enables  one  to  determine  whether  a  given  approximation, 
consisting  of  a  basic  scheme  together  with  corresponding  boundary  conditions, 
is  stable  or  not.  Roughly  speaking,  according  to  this  criterion  we  have  stability 
if  (and  only  if)  no  linear  combination  of  powers  of  roots  satisfying  a 
characteristic  equation  which  is  determined  by  the  basic  scheme,  may  serve  as  a 
non-trivial  solution  of  some  characteristic  boundary  constraints.  Thus,  in  order 
to  assure  stability  by  applying  the  above  stability  criteria,  one  must  first 
study  the  behavior  of  roots  of  the  corresponding  characteristic  equation.  This 
characteristic  equation  is  a  polynomial  equation  with  N  «  N  matrix  coefficients, 

N  denoting  the  order  of  the  original  approximated  system,  and  whose  degree 
depends  on  the  number  of  spatial  mesh  points  that  the  basic  scheme  rests  on. 
Studying  the  behavior  of  the  roots  of  such  an  equation  as  part  of  applying  the 
above  stability  criteria  for  general  difference  schemes,  is  a  complicated  task 
which  cannot  always  be  carried  out  fully  analytically.  Therefore,  examples 
given  in  the  literature  of  verifying  stability  for  initial-boundary  approximations, 
are  in  most  cases  restricted  to  specific  scalar  3-point  schemes. 

This  fact  motivates  us  to  look  for  simpler  sufficient  stability  tests. 
Scheme-independent  stability  tests  which  are  exclusively  dependent  on  the 
boundary  conditions,  are  particularly  useful  for  such  purpose.  Tests  of  this 
sort  have  two  main  advantages:  first,  their  being  independent  of  the  basic 
scheme  relieves  us  of  the  technical  difficulties  associated  with  the  computation 
of  roots  of  the  characteristic  equation  and  hence  the  procedure  of  checking 
stability  becomes  much  shorter;  and  secondly,  the  acquired  stability  is  not 


restricted  to  a  specific  approximation  but  instead  is  valid  for  a  family  of  basic 


schemes  which  are  characterized  by  some  genera.1  property. 

The  search  for  such  scheme-independent  stability  criteria  for  difference 
approximations  to  hyperbolic  initial -boundary  value  systems  is  the  main  subject 
of  the  dissertation. 

As  a  model  problem  for  the  general  linear  case  we  consider  the  hyperbolic 
system  u^  =  Au^  +  F  which  together  with  appropriate  initial  and  boundary 
conditions  is  well-posed  in  the  quarter  plane  x  >,  0,  t  >,  0.  In  prescribing 
these  appropriate  boundary  conditions  it  turns  out  that  one  must  distinguish 
between  inflow  and  outflow  (characteristic)  unknowns,  where  the  inflow  boundary 
values  have  to  be  determined  by  reflection  of  the  outflow  ones.  Based  on  this 
distinction  we  first  introduce  in  Chapter  1  a  general  method  of  numerical  boundary 
treatment  of  arbitrary  degree  of  accuracy,  such  that  the  entire  vector 
approximation  is  stable  if  and  only  if  the  scalar  components  of  its  outflow  part 
are;  thus  reducing  the  stability  question  to  that  of  a  scalar  (outflow)  problem. 
Therefore  from  that  point  on  our  discussion  concentrates  on  the  general  scalar 
approximation  as  it  is  represented  in  the  second  part  of  Chapter  1. 

In  Chapter  2  we  begin  the  stability  study,  drawing  on  the  stability  theory 
of  Gustafsson  et.  al.  [6]  which  we  briefly  survey  in  the  first  section  of  that 
chapter.  The  main  stability  criterion  in  that  theory  is  given  in  terms  of 
eigenvalues  and  generalized  eigenvalues  of  the  problem.  Then,  upon  reintroducing 
these  concepts  in  a  less  formal  manner  and  operating  under  the  four  basic 
assumptions  corresponding  to  those  which  were  made  in  [63,  we  may  apply  the  above 
criterion  [6,  Theorem  5-l3  which  states  that  a  given  initial -boundary  approximation 
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is  stable  if  and  only  if  it  has  neither  eigenvalues  nor  generalized  eigenvalues 
outside  the  unit  disc. 

In  the  second  section  of  Chapter  2  we  follow  the  analysis  in  [6]  which 
leads  to  the  formulation  of  the  main  stability  criterion  as  a  corresponding 
determinant  condition.  Then,  by  using  in  the  above  analysis,  a  suitable 
representation  of  the  general  form  of  (generalized)  eigenvalues  of  the  problem, 
we  derive  an  explicit  interpretation  of  the  determinantal  stability  criterion 
mentioned  above.  This  result,  which  seems  to  be  of  independent  interest,  is 
essential  for  the  general  stability  analysis  which  is  carried  out  afterwards 
and  is  needed  to  obtain  the  stability  criteria  of  the  desired  type. 

Chapter  3  —  the  main  one  in  this  work  —  discusses  scheme-independent 
stability  criteria.  The  boundary  conditions  considered  are  of  translatory 
type,  i.e.,  determined  at  all  points  in  the  boundary  domain  by  the  same 
procedure.  We  first  show  in  Section  3.1  that  when  dealing  with  such  boundary 
condition,  the  determinental  stability  criterion  obtained  in  Chapter  2  is 
equivalent  to  a  corresponding  scalar  condition.  This  scalar  condition  plays 
the  central  role  in  proving  the  scheme-independent  results  at  which  we  are  aiming. 

In  the  remainder  of  Chapter  3  we  3tate  our  main  results,  namely,  sufficient 
scheme-independent  stability  criteria.  These  results  are  obtained  upon  making 
two  quite  non-restrictive  assumptions  complementing  the  first  four  already  made 
in  Chapter  2.  It  is  shown  that  these  new  additional  assumptions  are  necessary  for 
our  scheme-independent  results  to  be  valid,  and  simple  scheme-independent  tests 
verifying  whether  a  given  problem  meets  these  assumptions,  are  provided.  We  study 
the  cases  of  both  one-level  and  multi-level  boundary  treatments.  In  the  (somewhat 
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simpler)  first  case,  the  well-known  result  (see  for  example  C 33 »  C T 3 )  stating 
that  two-level  stable  dissipative  schemes  together  with  (one-leveled) 
extrapolated  outflow  boundary  values  remain  as  conditionally  stable,  follows 
easily.  We  show  however,  that  this  widely  used  result  is  no  longer  valid  when 
dealing  with  multi-level  dissipative  schemes  involving  more  than  two  time-steps, 
unless  further  restrictions  are  made  on  them.  In  studying  the  wide  class  of 
multi-level  boundary  treatments,  we  will  employ  the  tools  of  dissipativity  and 
the  von  Neumann  condition  usually  used  only  in  connection  with  the  basic  scheme. 

The  stability  criteria  in  the  multi-level  case  are  given  in  terms  of  these 
concepts  which  are  well-understood  from  the  theory  of  pure  initial  value  problems 
and- whose  validity  can  be  easily  checked.  We  prove  that  an  arbitrary  stable 
dissipative  scheme  when  complemented  by  outflow  boundary  conditions  satisfying 
the  von  Neumann  condition ,  remains  stable  We  also  show  that  if  the  outflow 
boundary  conditions  are  dissipative ,  then  the  entire  approximation  is  stable 
independently  of  the  interior  scheme  (be  it  dissipative  or  a  non-dissipative  one). 

Finally,  in  Chapter  U,  we  utilize  the  above  scheme-independent  stability 
criteria  to  verify  the  stability  of  various  (outflow)  translatory  boundary 
conditions.  The  examples  considered  indicate  that  an  arbitrary  stable  dissipative 
scheme  whose  outflow  boundary  values  are  translatorily  computed  by  oblique 
extrapolation,  by  the  Box-Scheme  or  by  the  stable  weighted  Euler  scheme,  constitutes 
a  stable  approximation.  We  also  study  boundary  conditions  which  are  generated 
by  the  (right-sided)  implicit  and  stable  explicit  Euler  schemes.  Both  boundary 
treatments  are  found  to  be  unconditionally  stable  in  the  sense  that  when  augmenting 
arbitrary  stable  basic  schemes,  they  always  maintain  stability.  We  close  Chapter 
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U  by  considering  approximations  to  the  two  space  dimensional  problem 

u  =  au  +  bu  in  the  quarter  space  Xi0,ti0,  -“<y<®.  The  stability 
i-  x  y 

analysis  in  that  case  is  based  on  Fourier  transforming  with  respect  to  y 
.with  dual  variable  n),  thus  obtaining  a  one  space  dimensional  problem  of  the 
type  analyzed  in  previous  chapters  with  n-dependent  coefficients. 
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1.  THE  DIFFERENCE  APPROXIMATION  TO  THE  HYPERBOLIC  SYSTEM 


1.1.  The  reduction  to  the  scalar  problem 

Consider  a  first  order  hyperbolic  system  of  partial  differential  equations 


(1.1a) 


3u(x,t)/3t  *  A3u(x,t)/3x  +  F(x,t) 


in  the  quarter-plane  x  z  0,  t  >,  0,  with  initial  conditions 


(1.1b) 


u(x,0)  *  f(x)  ,  0  $  X  <  "  . 


Here,  u(x,t)  =  ( u^  ^ ^  (x,t ) , . . .  ,u^ (x,t ) ) '  is  the  transposed  vector  of  unknowns, 

A  is  N  x  N  non-singular  constant  coefficient  matrix  and 

F(x,t)  =  (F^(x,t),..,,F^(x,t))',  f(x)  i  (f^(x) . f(N)(x))'  ,  are 

N-dimensional  vector  functions. 

The  hyperbolicity  of  the  system  (1.1a)  implies  that  A  can  be  diagonalized 
by  a  similarity  transformation,  hence  we  may  assume  without  restriction  that  A 
is  already  given  in  its  diagonal  form 


(1.1c) 

I 


We  are  interested  in  the  uniqueness  of  the  solution  for  the  system  (1.1a). 


For  that  reason,  let  us  consider  the  partition 
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corresponding  to  that  of  A.  Since  the  characteristic  lines  associated  with  the 

N-L  components  of  u+(x,t)  go  outside  from  the  region  x,t  i  0,  it  follows  that 
+ 

u  (x,t)  which  is  carried  by  these  characteristics,  is  uniquely  determined  by  the 
initial  values  f(x)  in  the  whole  quarter-plane  x,t  *  0.  Because  of  the 
direction  of  its  characteristics,  *  u+(x,t)  is  considered  as  the  outgoing  part  of 
the  solution  u(x,t). 

The  L  characteristic  lines  associated  with  u  (x,t)  have  a  positive  slope 
and  hence  go  into  the  quarter-plane  x,t  >,  0.  Therefore,  in  order  that  u~(x,t), 
which  is  considered  as  the  ingoing  part  of  the  solution  y(x,t),  be  uniquely 
determined  in  the  quarter-plane,  it  is  necessary  to  specify  its  values  on  the 
boundary  line  x  =  0.  Thus,  for  the  solution  of  (l.la)  to  be  uniquely  determined, 

i 

we  prescribe  boundary  conditions  of  the  general  form 

(l.ld)  u“(0,t)  =  Su+( 0,t )  +  g(t)  ,  t  i  0, 


which  determine  the  missing  ingoing  boundary  values  by  reflection  of  the  outgoing 
ones.  Here,  S  is  an  L*(H-L)  constant  matrix  and  g(t)  =  (g^ (t ) , . . .  ,g  (t) ) ' 
is  an  L-dimensional  vector  function. 

To  solve  the  initial-boundary  value  problem  vi-.l)  By  a  difference  approximation 
ve  introduce  a  mesh-size  h  =  Ax  >  0,  At  >  0  such  that  X  =  At/Ax  ■  constant. 

Using  the  standard  notation  vy(t)  =  v(vh,t),  we  approximate  (l.la)  by  a 
consistent,  two-sided  multi -level*  scheme 


t 


(1.2a) 


Q_lvv(t+At)  = 


a=0 


Q  v  (t-aAt) 

<y  V 


At-Fv(t), 


=  1,2,. 


with  initial  values 


(1.2b) 


v^(aAt)  *  fy(aAt),  \>  >,  -r+1,  a  =  0,1, ...,s. 


Here 


Ev 

v 


v+1 


» 


are  difference  operators  with  matrix  coefficients  depending  on  A  and  on  A. 

In  order  to  determine  uniquely  the  solution  of  (1.2),  we  must  specify  at 
each  time  step  the  r  boundary  values  v^(t),  u  =  0,-1,.. ., -r+1.  For  the 
approximated  outgoing  unknowns  v*(t)  =  (v^+^(t),. . .  ,v^^(t))'  we  do  it  by 

boundary  conditions  of  the  form 


T 

S^v+(t+At)  *  /_  S^v+(t-oAt)  +  At*F  (t)  ,  u  *  0,-1,. . .  ,-r+l, 
-!  W  ^0  0  W  W 


(1.3) 


s<p)  »  5^CU)EJ  a  -  1  0  1  t  T  >  1 

Sc  CJa  E  *  -1,0,1,. .. ,r  ,  t  *  -1, 


where  are  (N-L)  *  (N-L)  constant  matrices  depending  on  A  and  on  A. 

J  ^ 

In  that  way,  we  maintain  the  property  Bhared  by  the  analytic  system  (1.1),  which 


is  the  unique  determination  of  the  outgoing  unknowns  in  the  whole  quarter-plane 
independently  of  the  ingoing  ones. 
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In  (1=3),  t  +  1  ( t  i  —1 )  denotes  the  number  of  time  levels  needed  for 
the  computation  of  v^(t+At).  The  case  t  =  -1  is  considered  to  be  the  case 
of  one-level  boundary  conditions  where  the  first  term  on  the  right  side  of  (1.3) 
is  taken  to  be  zero. 

For  the  computation  of  the  boundary  values  of  the  approximated  ingoing 
unknowns  v~(t)  =  (v|^  (L) , . . .  ,v^  (t) ) '  we  use  the  analytic  boundary  condition 

(1. Ua)  -  v“(t)  =  Sv*(t)  +  g(t). 


together  with  r-1  additional  conditions  of  the  form 


(1.4b) 


v  (t) 

w 


-t » 

j=_r+l 


UJvJ(t)  .  eu(t) 


U  *  -1,-2,. . .  ,-r+l. 


Here,  are  Lx(N-L)  constant  matrices  and  the  g^(t)  are  L-dimensional 

bounded  vector  functions  depending  on  h  and  on  g(t).  In  other  words,  as  in 
the  analytic  case,  the  computation  of  the  ingoing  boundary  values  is  based  on 
reflection  of  the  outgoing  ones. 

It  is  well-known  that  using  conditions  of  the  general  form  (1.3),  one  can 
achieve  at  the  boundary,  arbitrary  degrees  of  accuracy.  We  note  that  this  is 
true  also  for  conditions  of  the  type  (1.4b).  Indeed,  if  accuracy  of  order  d 
is  desired  we  can  use  the  Taylor  expansion  of  a  smooth  solution  for  (l.l) 


u^(t) 


J=0  3‘  dxJ 


I 
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and  by  the  differential  system  (1.1a)  and  (l.ld)  we  get  for  a  typical  spatial 
derivative  in  the  above  expansion 

t)  =  (A")  J  -^u-(0,t)  =  (A“)~J[S-4u+(0,t)  +  -^«(t)]  = 

3xJ  3tJ  3tJ  dtJ 

=  (A')_J[S(A+)J— 4  u+(0,t)  + 

3xJ  dtJ 

Thus,  (1.4b)  follows  upon  approximating  3^/3x^u+( 0,t )  by  linear  combinations 

of  u+  (t) , . . . ,u+(t)  of  the  right  accuracy. 

-r+1  q. 

The  difference  approximation  is  completely  defined  now  by  the  (basic)  scheme 
(1.2a)  together  with  the  boundary  conditions  (1.3),  (1.4)  and  we  raise  the 
question  of  its  overall  stability  which  means,  according  to  Definition  3-3  in 
[6],  that  the  discrete  solution  v^(t)  could  be  estimated  with  the  aid  of  the 
inhomogeneous  terms  F^(t) ,  v  }  -r+1,  g(t)  and  g^(t),  M  *  -1,-2,. . . ,-r+l. 

For  that  purpose,  we  split  the  scheme  (1.2a)  into  its  inflow  and  outflow 

parts 

(1.5)  Q^v'U+At)  =  'y  Q~v~(t-aAt )  +  At*F~(t),  v  =  1,2,..., 

o®0 

s 

(1.6)  Q^v^t+At)  =  'y  Q*v*(t-cAt)  +  At*F*(t),  v  *  1,2,..., 

a*0 

which  are  coupled  through  the  boundary  conditions  (1.4).  Here,  Q~,  Q+  are 

a  o 

difference  operators  which  are  given  by 


i 

i 

> 


I 
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Thus,  Q^,  Q*  denote  respectively  the  partition  of  the  difference  operators 

into  their  inflow  and  outflow  parts,  according  to  the  dependence  of  the  matrix 

—  + 

coefficients  A,  on  A  and  A  . 

Jo 

In  order  to  assure  the  stability  of  the  entire  approximation,  both  of  its 

parts  -  the  outflow  part  (1.6),  (1.3)  and  the  inflow  part  (1.5)  ,  (l.M  -  have  to 

be  stable.  We  note  that  the  outflow  approximation  (1.6),  (1.3)  is  independent  of 

the  inflow  values,  while  the  inflow  approximation  (1.5) »  (1.1*)  depends  on  the 

outflow  part  only  to  the  extent  that  the  outflow  computations  provide  the 

inhomogeneous  boundary  values  in  (l.Ub). 

Let  us  consider  first  the  stability  of  the  self-contained  outflow  approximation 

(1.6),  (1.3).  Since  the  difference  operators  Q*,  are  expressed  in  terms 

of  A*  and  c[y  ^  which  in  turn  depend  on  the  diagonal  matrix  A+ ,  it  follows 
Ja  Jo 

that  the  outflow  problem  splits  into  N-L  independent  scalar  approximations , 

Thus,  -t-he  outflow  problem  is  stable  if  and  only  if  its  N-L  scalar  component 
approximations  are. 

Now,  suppose  the  outflow  approximation  was  found  stable,  then  it  remains  to 
determine  whether  the  same  holds  true  for  the  inflow  approximation  (1.5),  (l.M. 
Since  the  out f lor. •  values  determined  by  the  stable  outflow  computation  are  bounded, 

then  the  summation  ^  D  .v*(t),  which  appears  on  the  right  side  of  (lJ‘b)  is  a 
J=-r+l  UJ  J 

bounded  term,  independent  of  the  inflow  computation.  Thus,  the  right  side  of 
(l.Ub)  consists  of  two  bounded  terms,  which  are  independent  of  the  inflow 
computation  and  therefore,  for  the  purpose  of  determining  stability,  it  may  be 
considered  as  an  arbitrary  inhomogeneous  term  that  provides  the  ingoing  boundary 


*> 


values . 


Recalling  that  the  difference  operators  Q-  defining  the  basic  inflow 

o 

scheme  are  expressed  in  terms  of  A.  ,  which  in  turn  depend  on  the  diagonal 

J  ® 

matrix  A  ,  it  follows  that  the  inflow  problem  splits  into  L  independent  scalar 
initial-boundary  approximations  the  boundary  values  of  which  are  determined  by 
some  inhomogeneous  bounded  terms.  Thus,  the  inflow  problem  is  stable  if  and 
only  if  its  L  corresponding  scalar  components  are. 

Concerning  the  stability  question  of  the  initial-boundary  approximation 
which  is  discussed  above,  it  is  obviously  necessary  to  require  the  stability  of 
its  basic  scheme  should  it  be  applied  to  the  pure  initial -value  problem, 

—  00  <  <  oo 

Assume  that  the  basic  inflow  scheme  indeed  satisfies  the  above  necessary 
stability  requirement.  Then,  as  we  shall  see  later  on,  it  follows  that  each 
scalar  component  of  the  inflow  approximation  whose  boundary  values  are 
determined  by  an  arbitrary  bounded  term,  is  unconditionally  stable,  [7],  ClO], 
and  we  therefore  obtain  the  unconditional  stability  of  the  entire  inflow  part, 

(1.5).  d.u). 

We  conclude  that  the  entire  approximation  is  stable  if  and  only  if  its 
scalar  components  are.  Furthermore,  according  to  the  remark  above,  it  is 
sufficient  to  consider  only  the  outflow  ones.  Thus,  in  both  cases  -  either 
the  inflow  case  or  particularly  the  outflow  one,  it  is  the  scalar  approximation 
the  stability  of  which  we  have  to  look  for,  so  hereafter  we  may  restrict  our 
discussion  to  the  scalar  approximation,  bearing  in  mind  that  our  forthcoming 


results  go  over  to  the  general  vector  case. 
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1.2  The  scalar  approximation  and  its  solvability 

Consider  the  scalar  hyperbolic  initial -boundary  value  problem 


(1.7a)  3u(x,t)/3t  =  a3u(x,t)/3x;  a  =  constant  4  0;  u(x,0)  =  f(x);  x,t  z  0. 


Whereas  the  outflow  problem,  a  >  0,  is  well-posed  in  L2(0,»),  the  inflow  problem, 
a  <  0,  is  not,  unless  suitable  boundary  conditions  are  given  at  x  =  0.  Therefore, 
we  shall  examine  (1.7a)  together  with 

-  no  boundary  conditions  for  the  outflow  problem,  a  >  0; 

(1.7b) 

-  boundary  conditions  u(0,t)  =  g(t),  t  }  0  for  the  inflow  problem,  a  <  0. 

To  approximate  (1.7)  numerically,  we  set  a  time  step  At  >  0  and  a  mesh 

width  h  =  Ax  >  0,  a  grid  function  v'('t)  =  v(vh,t),  v  =  0,  ±1,  +2,...,  and  a 
consistent  multi-level  finite  difference  scheme 


(1.8) 


a 

Q_1vy(t+At)  =  y  Q^v  (t-oAt);  v  *  l,2,...,t  z  sAt 
o=0  0  V 

% =  %r  VJ  !  0  ■  . . .  e3vv  "  Vr 


Here,  r,p  >  0  and  s  are  natural  numbers  and  the 
depend  on  the  coefficient  a  and  the  fixed  ratio 


a.  's  are  constants  which 

Jo 

X  r  At/Ax  ■  constant. 


We  note  that  the  consistency  of  scheme  (1.8),  l.e.,  its  being  at  least 


first  order  accurate,  may  be  characterized  by  the  following  two  equations: 
zero  order  accuracy  requirement 


[1.9a) 


J=-r  A  aa0  J*-r 


and  the  additional  requirement  for  first  order  accuracy 


(1.9b) 


s  ^£_  a 

/  Ja.  ,  =  2_  Z_  Ja1fI  -  la  2__  (°+l)  / 


a*0  J=-r 


The  equalities  (1.9a),  (1.9b)  may  be  written  respectively,  in  the  following 


compact  form 


(1.10a) 


(1.10b) 


i  j,ju)iz.i  ■  -x*  ;  [  v  ’  ^  • 


where  the  scalar  functions  a,(z)  are  defined  by 

V 


(1.11) 


8 

‘jU)  •S*'*'1*!..  *  *3.-1  •  'r 


$  3  $  P* 


and  following  [63,  [83,  we  shall  operate  under 
ASSUMPTIOH  I  (Assumption  5-5  in  [63). 


a  „(z),  a  (z)  *  0  ,  |z|  i  1. 

««F  P 


* 
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It  is  clear  that  under  Assumption  I,  the  vector  coefficients. 

(a  , ,a  a  )'  does  not  vanish,  and  since  r  >  0  it  follows  that  in 

— 3?  $  ”1  “1*^0  **r  y  s 

order  to  assure  the  uniqueness  of  the  solution  of  (1.8),  we  have  to  supply  its 
discrete  values  at  the  boundary  points  x^,  p  =  0,  -1, ...,-r+l.  These  will  be 
defined  via  boundary  conditions  of  the  form 


(1.12a) 


T 

S^vy(t+At)  *  y~  S^v^(t-aAt);  V  =  0,-l,.~,-r+l,  t  5  xAt, 

An  E>)  i  0  =  -1*0 . .  T  *  -1- 

0  J^O  J° 


Here,  the  c|U^'s  are  constants  which  depend  on  a  and  X,  q  is  a  natural 
number  and  t+1,  t  i  -1,  indicates  the  number  of  previous  time  levels  which  we 
need  in  order  to  compute  the  boundary  values  at  the  next  time  level,  t+At.  We 
note  that  the  one-level  boundary  conditions  namely  t  =  -1,  is  a  special  case  of 
(1.12a)  whereupon  (1.12a)  takes  the  form 


(1.12b) 


(t+At)  =  0, 

-1  p 


0,-1,. . . ,-r+l,  t  i  0. 


It  is  clear  that  the  computation  of  the  boundary  values  via  the  r  boundary 
relations  (1.12),  the  linear  independence  of  which  is  assured  by  taking 


(1.13) 


ci,*^..  0  ,  p  =  0,-1,. ..  ,-r+l. 


is  done  in  the  specified  order,  namely,  p  *  0,-1, ... ,-r+l. 


Now,  the  basic  scheme  (1.8)  together  with  the  boundary  conditions  (1.12) 


I 


-In¬ 


completely  define  the  finite  difference  approximation,  whose  numerical  solution 
is  initiated  with  the  aid  of  the  initial  values  given  by 


( 1.14 ) 


v^(oAt)  =  f^(adt)  ,  v  >,  -r+1,  a  =  0,1,..., s. 


Following  [6],  we  define  the  solvability  of  the  difference  approximation  as 
the  property  of  being  able  to  uniquely  obtain  bounded  grid  values  at  t+ht  by 
applying  (1.8)  and  (1.12),  thus  making  use  of  the  discrete  values  which  were 
already  computed  at  previous  time  levels. 

To  ensure  solvability,  we  consider  the  space  l^x) ,  of  all  the  grid 


functions  w 


(w  i 

v  v=-r+l 


satisfying 


00 


Z_  (wvf2  <  -• 

vs-r+i 


Upon  defining  respectively 


an  inner  product  and  a  norm  by 


(v,w)x 


Ax •  /  v  w 

‘  v  v 
v=-r+l 


(v,w)x. 


l^(x)  becomes  a  Hilbert  space,  a  discrete  analogue  to  L2(0,®).  Now,  denote  by 
w  .  v  i  -r+1,  the  discrete  values  to  be  computed  at  the  next  time  level,  t+dt 

v  » 

and  rewrite  the  approximation  (1.8),  (1.12)  in  the  form 

(1.15)  Q  ^wy  •  v  ■  1,2,...;  s^wjj  =  ’•'y  V  *  0,-1,. . .  ,-r+l, 

where  ?■{?)"  .  T  €  t-(x),  stands  for  the  linear  combinations  of  previously 

v  v=-r+l  2 
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computed  values,  as  given  on  the  right  sides  of  (1.8)  and  (1.12). 

The  solvability  of  the  approximation,  which  will  be  henceforth  assumed 
throughout  this  work,  is  thus  cast  in  the  following  form. 

ASSUMPTION  II  (Assumption  3.1  in  [6]  ). 

There  exists  a  constant  KQ  >  0,  such  that  for  every  V  £  l^ix)  there  is  a 
unique  solution  w  £  ^(x)  for  (1.15)  with 


2 

x‘ 


We  note  that  the  solvability  condition  is  automatically  fulfilled  in  case 
that  the  basic  scheme  (1.8)  is  explicit,  i.e.,  Q  ^  e  constant* I.  Concerning 
the  solvability  m  the  general  implicit  case,  the  following  result  due  to  Osher 
[12]  (based  on  Strang's  earlier  paper  [l6])  holds. 

LEMMA  1.1. 

Let  the  index  rQ,  0*rQsr  be  defined  by 


(1. 16 ) 


rQ  *  max{j|a_j  _1»‘0,0$J«r] 


and  let  J  *  1,2,. . . ,m  be  the  m  roots  counted  according  to  their 

J 

multiplicities,  of 

(1.17)  Q  -.(*)  ■  '/Li  a,  , le*5  =  0, 


which  are  lying  inside  the  unit  disc  0  <  |<j|  <  1. 
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I •  ( Theorem  I  in  [ 12 ] ) .  The  following  three  conditions  are  necessary  and 
sufficient  for  solvability: 


(1. l8a) 


Q_l(e  >  *  0  ’  l?l  5  w 


( 1 . 18b ) 


(l.l8c)  the  associated  problem  with  (1.15)  which  consists  of  the  basic  scheme 


(1.19a) 


i  a1  _i 


Vj  =  0,  v  =  1,2 . 


together  with  the  ii 


aoundary  conditions 


(1.19b) 


S  M/w  =  V  ,  u  =  0,-1 . -r+1, 

-1  li  w 


has  a  unique  solution  in  £2(x). 


II.  Let  the  basic  scheme  (1.19a)  be_  a  right-sided  onc,i.e. ,  r^  *  0.  Then 


(1.20) 


Q_x(<)  4  0  ,  0  <  M  s  1 


is  a  sufficient  condition  for  solvability. 

In  particular  solvability  follows  for  explicit  basic  schemes,  where 
we  have  rQ  a  0,  Q  ^(k)  =  constant. 


I 


f 
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PROOF ,  By  (l,l8a)  we  may  apply  the  Argument  Principle  for  Q  ^(k)  on  the 

unit  circle  |k|  =  1,  obtaining  that  (l.lffb)  is  valid  if  and  only  if 
ir 

(2n)  d[argQ  .(e^)]  =  0.  Hence,  conditions  (1.18)  are  exactly  those  of 

-IT 

Osher  ((d)  (e)  and  (g)  in  [12])  and  by  Theorem  I  in  [12]  they  are  equivalent 
to  solvability. 

To  prove  the  second  part  of  the  lemma  we  first  note  that  (1.20)  implies 
(l.l8a).  Also  (1.20)  implies  that  Q  ^(k)  =  0  has  no  solutions  inside  the 
unit  disc,  i.e.,  m  =  0,  and  since  r^  =  0  we  have  (l.l8b)  as  well.  Finally 
since  by  (1.20)  Q  1(x)  =  0  has  no  solutions  in  the  closed  unit  disc,  it 
follows  that  the  most  general  solution  of  (1.19a)  in  *2(x)  must  vanish,  i.e., 
w^  =  0,  v  =  1,2, . . .  .  In  addition  the  boundary  values  w^,  u  *  0  ,-l , . . . ,-r+l , 
can  be  uniquely  computed  by  applying  ( 1.19b)  in  the  successive  order 
w  =  0,-1, ... ,-r+l ,  so  we  get  (l.l8c).  Having  (l.l8a),  (l.l8b),  (l.l8c),  part 
I  of  the  lemma  completes  the  proof. 


I 
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2.  STABILITY  ANALYSIS 


f? . 1 .  The  stability  definition  and  Gustafsson's  et .  al.  Main  Theorem  C6] 

In  a  similar  way  to  the  above  definition  of  the  space  l^(x) ,  we  introduce 
the  discrete  spaces  £  (t)  and  l^(x,t) ,  which  become  Hilbert  spaces  upon 
defining  respectively  an  inner  product  and  a  norm  by 


(v 


(v,w)  =  At-y1  v(oAt)w(oAt)  ;  p f  =  (w,v)  , 

t  a=0  t  t 

OO  CD 

,w)  =  At*Ax*2  / _  v  (oAt)w  (oAt);  !!w!2  =  (w,w)  . 

X*t  £o^r+lv  v  ’  x,t  x,t 


Now,  let  us  write  the  difference  approximation  (1.8),  (1.12),  in  the 

operational  form  ,  > 

max\ t , s ) 

(2.1 )  G  ^v( t+At )  =  G  v(t-oAt);  v(t-oAt)6t2(x) , 

a=0  0 


where  — >  l^(x)  are  linear  bounded  operators  determined  by  the  basic 

scheme  (1.8)  together  with  the  boundary  conditions  (1.12).  Here,  the  solvability 
assumption  II,  is  expressed  by  the  fact  that  G  ^  has  a  bounded  inverse  in  the 
whole  of  fcg(x). 

DEFINITION  2.1  (Definition  3.3  in  [6]). 

Consider  the  inhomogeneous  approximation  associated  with  (2.1) 

max(T.s) 

(2.2)  G_1w(t+At)  -  2  0  w(t-oit)  ■  At*F  F  =  £  igU), 

o*0 

together  with  vanishing  initial  values  f^(oAt)  ■  0.  The  approximation  (1.8), 


TT> 


/ 
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(1.12)  is  said  to  be  stable  if  there  exist  constants  >  0  and  >,  0, 

oo 

such  that  for  every  F and  every  a,  a  z  a^,  the  solution  w  =  r+^» 

of  (2.2)  satisfies  the  estimate 


(2.3) 


< 


2 


Pr>^  b 


At‘e 


-a(t+At] 


U=-r+l 


f  !2 

u't 


♦  I, 


t-a(t+At'  {2 


x  ,t 


It  is  of  course  understood  that  before  turning  to  investigate  the  stability 
of  the  initial -boundary  approximation  (2.1),  one  has  to  assure  first  the  stability 
of  scheme  (1.8)  should  it  be  applied  to  the  pure  initial-value  problem,  -°°<x<”. 
ASSUMPTION  III  (Assumption  5-1  in  [6]). 

The  scheme  (1.8)  is  a  stable  approximation  for  the  Cauchy  problem, 

It  is  well-known  (see  for  example  [13,  Section  U])  that  Assumption  III, 
may  be  characterized  by  the  two  following  conditions: 

(i)  The  von  Neumann  condition;  namely,  the  z-solutions  for  the  eigenvalue 
problem 

(2. ba)  Q  ,(U)  -  z"0"1^  (1C)  =  0,  QUO  3 

a=0 


which  may  be  rewritten  in  the  equivalent  form 


(2.4b) 


a .(z)e1^  ■  0;  a,(z)  =  -5~~ z~°~ 

J»-r  i  3  o«0 


aJa  +  aJ,-l’ 
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satisfy  |z  =  z(C)|  i  1  for  |sj  S  * 

(ii)  The  solutions  for  the  eigenvalue  problem  (2.1*)  which  are  lying  on  the 
unit  circle,  are  simple,  i.e.,  for  z^  =  Zq(5q),  |zq|  =  1,  satisfying  (2.M, 
we  have 


(2.5) 


/Li  a1'(z)e 

J=-r  J 


iJCr 


i„_  i  o, 

z"zo 


[  ]• 


d[  ] 
dz 


In  addition  we  require 

ASSUMPTION  IV  (Assumption  5.4  in  [6]). 

Denote  by  z,(0  the  solutions  fcr  (2.U).  Then,  the  scheme  (1.8)  is 
J 

either  dissipative,  i.e., 

bj(£)|  <  1  .  0  <  |t|  J  n, 


or,  it  is  nondissipative,  i.e.. 


|ZjU)|  =1  ,  0  s  |F,|  S  it. 

Operating  with  our  scalar  approximation  under  Assumptions  I-IV,  enables  us 
to  use  the  results  obtained  by  Gustafsson  et,  al. [6]j  in  particular,  we  are 
interested  in  their  main  result,  characterizing  the  stability  of  an  initial¬ 
boundary  approximation.  The  remainder  of  this  section  is  therefore  devoted  to 
a  brief  survey  of  some  of  the  points  concerning  this  matter. 

Consider  the  z-eigenvalue  problem  ~  the  associated  resolvent  equation  — 

6iVen  by  max(x.s) 

(2.6)  G(z)*  «  0,  G(z)  =  G  -  > 

o*0 
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/ 


which  follows  upon  substituting  a  grid  solution  of  the  form  v^(t)  =  2^^$^ ,vi-r+l , 
into  the  approximation  (2,1). 

DEFINITION  2,2, 

A  complex  number  Zq,  is  said  to  be  a  spectrum-point  of  the  approximation 
with  associated  eigenvector  *,  $  ?  0,  if  there  exists  a  sequence  of  vectbrs 

€  l2(x),  satisfying 


(2.7) 


G(zq)* 


(3)  I  I 


0,  * 


(3)  JLL 


♦  ;  J  — > 


We  note  that  the  eigenvector  ♦  ,  associated  with  the  spectrum  point  zQ,  is 
not  necessarily  in  ^^(x).  In  case  the  vector  $  is  indeed  in  the  point 

zQ  is  an  eigenvalue  of  the  approximation;  otherwise,  when  *  ?  l2(x),  zQ  is  a 
generalized  eigenvalue  of  the  approximation.  In  either  case,  the  boundedness  of 
the  operator  G(zQ)  implies  that  the  corresponding  eigenvector  *  satisfies 

(2.8)  G(zQ)*  =  0. 

It  is  not  hard  to  see  that  in  order  to  assure  the  stability  of  the  difference 
approximation,  we  have  to  assure  first  that  the  necessary  condition  of  Ryabenkii- 
Goudunov  is  to  be  fulfilled,  namely,  that  the  approximation  has  no  eigenvalues  z 
with  |z|  >1  (Lemma  4.1  in  [6]}.  Indeed,  the  existence  of  such  an  eigenvalue 
Zq,  I*  |  >  1,  with  associated  eigenvector  ♦,  ♦  €  l^ix) t  implies  that  the  grid 
function 

(2.9)  vv(t)  -  zJ/At  fv  ,  v>,  -r+1 , 


I 
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is  a  solution  for  the  approximation  (1=8),  (1.12),  with  initial  values 
v^(aflt)  =  v  2  -r+1,  o  =  0,1,..., s,  a  solution  which  exponentially  diverges 

with  the  refinement  of  the  grid  as  At  — >  0.  Evidently,  such  a  divergence 
cannot  he  allowed  within  the  limits  of  any  stability  definition  and  in  particular 
definition  2.1,  (Theorem  3.1  in  [6]). 

The  main  result  in  [6]  strengthens  the  necessary  condition  mentioned 
above  to  be  also  sufficient. 

THEOREM  2-1.  (Theorem  5=1  in  [6]). 

The  difference  approximation  (1.8),  (1.12)  is^  stable ,  if  and  only  if  it  has 
neither  eigenvalues  nor  generalized  eigenvalues  z,  with  | z|  i  1. 
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2.2.  A  determinantal  stability  condition 

In  this  section,  we  intend  to  express  the  stability  condition  which  was  given 
in  Theorem  21,  in  a  suitable  algebraic  formulation. 

For  that  purpose,  we  consider  the  characteristic  equation ,  associated  with  the 
basic  scheme  (1.8), 


(2.10) 


P(z 


»<)  =  Z__  a, 

j  =  -r  J 


(z)«cJ  =  0, 


whose  r+p  roots  counted  according  to  their  multiplicities,  are  continuous 

functions  of  z. 

The  behavior  of  these  z-dependent  solutions ,  plays  a  central  role  in 
determining  the  set  of  the  spectrum  points  of  approximation  (1.8),  (1.12). 

The  following  lemma  summarizes  the  results  which  were  given  in  [6  ]  and  [8  ] 
concerning  those  solutions,  for  any  solvable  approximation. 

LEMMA  2.1. 

Consider  the  solvable  approximation  (1.8),  (1.12). 

I.  (Lemma  5.2  in  [6]).  When  scheme  (1.8)  satisfies  the  von  Neumann  condition, 

then  the  r+P  solutions  of  its  associated  characteristic  equation  (2.10) 
are  split  for  |z|  >  1: 

£  with  |  ( z )  |  >  1,  and  the  rest  r  solutions  with  0  <  |ic^(z)|  <  1. 

II.  (Lemma  2  in  [8]).  When  Bcheme  (1.8)  is  of  dissipative  type  and  it 

additionally  satisfies 

(2.11)  S  a  (z-eiiP)  4  0  ,  0  <  |<p|  $  »  , 

J— r  J 

then  the  above  splitting  property  is  maintained  for  )z|  s  1,  z  I*  1. 


I 
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PROOF  The  proof  of  both  parts  of  the  theorem  is  based  on  the  idea  of  identifying 
the  z-values  for  which  the  characteristic  equation  (2.10)  has  no  ic-solutions  on 
the  unit  circle,  i.e,,  tc  =  e1^,  0  $  |c|  s  if. 

We  first  note,  that  the  solutions  <^(z),  |z|  >„  1  of  (2.10),  are  exactly 

those  which  solve 


(2.12) 


<r*P(z,tc) 


a,  ( z)ie^+r  =  0; 
J— r  J 


indeed,  multiplying  the  characteristic  equation  (2.10)  by  the  factor  <T  as  in 
(2.12),  does  not  yield  additional  zero  solutions,  <  -  0,  since  by 
Assumption  I, 

<r«P(z,K)|(c=0  *  a_r(z>  *  °»  M  *  1* 


Now,  concerning  the  first  case,  the  von  Neumann  condition  implies  that  the 
characteristic  equation  (2.12)  has  no  K-solutions  on  the  unit  circle  for  all 
z-values  with  |z|  >  1,  since  by  (2.4b)  we  have  for  |z|  >1 

(2,13)  eir^*P(  z.ic^e1^)  =  "  a,(z)ei^+r^  i  0,  0  f  |  C I  S 

JB-r  J 

Concerning  the  second  dissipative  case,  a  slight  change  is  needed  in  Kreiss' 
original  proof  [8,  Lemma  2]  in  order  that  the  splitting  property  will  be  valid 
for  our  mu lit-level  scheme  (1.8). 

In  this  case,  the  dissipativity  property,  which  indicates  that  the  solutions 
zU)szU«eU)  of  (2.12)  satisfy 


(2.14) 


|z(C)|  <1  ,0«  |C|  *i, 


•  r 


■tr 


I 
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implies  that  for  |z|  i  1,  we  have 


(2. 15) 


r?-P(2,K=eU)  s  YL  *  0  ,  0  <  U|  S 


j=rJ 


It  therefore  remains  to  check  the  single  point  £  *  0. 

By  continuity  considerations,  it  follows  from  (2.lU)  that 


(2.16) 


|zU=0)|  $  1. 


We  recall  our  hypothesis  (2.11),  which  implies  that  the  solutions 
z  =  z(£)  of  (2.12)  for  ic  =  el5|?=0  satisfy 

(2.17)  zU=0)  4  e1(P  ,  0  <  M  5  *  . 

Thus,  by  (2.l6),  (2.17),  we  obtain  for  |z|  *  1,  z  4  1 

(2.16)  eirC*P(z,K»eU)|^0  =  ^  «jU>  4  0. 

Combining  (2.15)  and  (2.18)  yields  that  (2.13)  is  valid  in  the  second 
case,  for  all  the  z-values  satisfying  |z|  >1,  z  4  1. 

Now,  the  solutions,  k±{z),  for  (2.12),  are  continuous  functions  of  z. 
Therefore,  the  number  of  solutions  ic^  satisfying  0  <  |k^(z)|  <  1,  is 
independent  of  z,  as  long  as  |z|  >  1  in  the  first  case,  and  |z|  *  1,  z  4  1, 
in  the  second  one.  By  letting  z  tend  to  infinity,  |z|  —>  •  ,  it  follows 


I 
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that  this  number  is  equal  to  the  number  of  solutions  tc  inside  the  unit  disc  of 


(2.19) 


cr*P(  |z|  — >  *,*)  =  icr- Q  .  (<)  =  'y  a  x^+r  = 


In  order  to  find  how  the  solutions  of  (2.19)  are  split  we  denote,  as  in 

(l.l6),  by  rn,  0  s  rn  5  r,  the  maximal  index  for  which  a  0,  and 

u  u  -ro>~1 

rewrite  (2.19)  in  the  form 


*  0. 


Then,  the  number  of  solutions  k,  0  $  |ic|  <  1  of  (2.19),  consists  of  r-rQ 
zero  solutions,  <  -  0,  and,  by  Lemma  1.1  which  is  valid  under  our  solvability 
assumption,  Tq  additional  solutions  of 


Hence,  there  exist  r  solutions  inside  the  unit  disc,  and  p  outside  it,  and 
the  result  follows. 

REMARK  2.1.  We  note  that  in  the  course  of  proving  Lemma  2.1,  we  didn't  need 
Assumption  I,  except  to  assure  that  multiplying  the  characteristic  equation 
(2.10)  by  the  factor  xr  ,  as  given  in  (2.12),  does  not  yield  additional  zero 
solutions,  k  ■  0. 

REMARK  2.2.  To  assure  the  splitting  property  for  |z|  *  1,  z  +  1  in  the 
dissipative  case,  then  according  to  the  second  part  of  the  above  lemma, 
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condition  (2  31)  has  to  hold.  We  note  that  this  condition  is  actually  part  of 
the  splitting  property  since  violation  of  (2.11)  means  that  for  some 
ipQ,  0  <  |tp0!  s  ir,  we  have 

P(z=eitf0,K=l)  =  'y  a  (z=eiy3°)  =  0; 

J*-r  J 


that  is,  if  (2.11)  is  violated  then  for  some  z  with  |z|  =  1,  z  t  1, 
the  characteristic  equation  (2.10)  ha3  a  root  (ic=l)  on  the  unit  circle. 
Moreover,  condition  (2.1l)  is  generally  necessary  in  the  sense  that  it  is 
independent  of  dissipativity.  Indeed,  let 


v^(t+At)  = 


t 

J=-r 


aJ,0Vv(t) 


v  =  1,2,.. 


I 


be  any  two-level  dissipative  scheme  which  is  (at  least)  zero  order  accurate, 
i.e.,  by  (1.10a)  we  have 


(2.21) 


y  a.(z»l)  ■  1 
J--r  J 


0. 


Now,  let  s  be  a  positive  integer  and  consider  the  solvable  scheme 


(2.22) 


v  (t+At) 
v 


v 


(t-sAt) 


v  *  1,2, 


I 


This  scheme  is  dissipative  yet  (2.11) 

.  2viJ/(s+l)  ,  ,  „ 

for  z  ■  Uj ,  Uj  *  e  ,  J  ■  1,2 


is  violated  since  by  (2.21)  we  have 
,  • .  •  ,s 


r 
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(2.23) 


2=0) 


J 


1—2 


-( S+l) 


J=-r 


J.O 


0. 


The  last  example  (2.22),  which  shows  the  necessity  of  the  additional  condition 
(2.11),  is  of  course  a  degenerated  one.  When  we  turn  to  examine  either  the  two- 
level  or  the  three-level  schemes,  i.e.,  s  =  0  or  s  =  1,  which  are  apparently 
the  ones  used  most  often,  we  find  that  the  additional  condition  (2.11),  may  be 
omitted,  or,  at  least  may  be  weakened.  This  is  the  content  of  the  next  lemma. 

LEMMA  2.2. 

Let  the  scheme  (1.8)  be  accurate  of  (at  least)  order  zero  (that  is ,  even 
the  consistency,  (1.10),  is  not  necessarily  required). 

I .  For  two-level  scheme,  s  ■  0 ,  we  have 


(2.24) 


i  a  (x=ellp)  +  0 
J=-r  J 


0  <  |<p|  S  Tt. 


II.  For  three-level  scheme ,  s  =  1,  which  satisfies 


(2.25) 


(z*-l) 


*  0, 


(2.24)  still  holds. 

Thus,  the  additional  condition  (2.11)  is  automatically  fulfilled  in  the 
two-level  case,  and  has  to  be  verified  at  the  single  point  z  ■  -1  in  the  three- 
level  one. 

PROOF.  When  scheme  (1.8)  is  of  (at  least)  zero  order  accuracy,  (1.10a)  implies 


•wc- 


t 


33  - 


(2-26) 


P(z=1,k=1)  =  2_  a. 

J-r  J 


(z=l)  =  0. 


Now,  in  the  two-leveled  case,  the  characteristic  function 


(2.27) 


P(z,ie=l)  =  Q_x(ic=l)  -  z-1'Q0(ic=l) , 


-1 


is  a  polynomial  of  lirst  degree  in  the  argument  z  ,  and,  by  (2.26),  its  only 


root  is  z  =1;  hence 


P(z=e 


i<p,ic=l)  =  ^  a  (z=elcf))  ¥  0  ,  0  <  |<p|  5  *, 


j=-r 


thus,  (2.24)  holds. 

In  the  three-leveled  case,  the  characteristic  function 


(2.28)  P(z,k=1)  =  Q_^(k*1)-z-^'Qq(»c=1)  -  z^Q-^ic*!) 


is  a  polynomial  of  second  degree  in  the  argument  z_\  whose  coefficients  are 
real,  and  by  (2.26),  z  1  »  1,  is  one  of  its  two  roots.  Hence,  the  other  root 
of  (2.28)  is  real,  and  therefore 


(2.29) 


P(  z*e^,K=l) 


=  a  (z=eil<>)  ¥  0  ,  0  <  |<p|  <  n 
J«-r  J 


How,  combining  (2.29)  with  our  hypothesis  (2.25)  which  merely  asserts 
that  (2.29)  is  valid  also  for  |<p|  ■  it,  gives  the  desired  result. 


r 
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We  return  now  to  the  characteristic  equation  (2.10),  the  solutions  of 
which  were  discussed  in  the  Lemma  2.1. 

Denote  by  =  k^(z),  the  distinct,  z-dependent  solutions  of  the  characteristic 
equation  (2.10)',  each  with  corresponding  multiplicity  m  =  m^z).  Since  our 
basic  scheme  (1,8)  is  always  assumed  to  satisfy  the  von  Neumann  condition 
(Assumption  III),  then  by  the  splitting  property  stated  in  Lemma  2.1,  we  may 
distinguish  between  two  groups  of  solutions  of  (2.10): 

the  group  of  the  inner  solutions,  k^(z),  lsfcsn,  which  are  characterized  by 

0  <  |iet(z)|  <1  ,  |z|  >  1; 

and  the  group  of  the  outer  solutions ,  containing  the  rest  of  the  solutions,  and 
characterized  by 

I^U)!  >  1  ,  |z|  >  1. 

Note  that  by  continuity ,  the  inner  (outer)  solutions  are  well-defined  for  \z\  >,  1, 

where  the  milder  inequalities  |k^|  $  1  4  are  valid. 

Now,  let  z  with  |z|  z  1  be  given.  If  z  is  an  eigenvalue  or  a 

generalized  eigenvalue  of  the  approximation,  then  there  exists  a  corresponding 

(nontrivial)  eigenvector  *  such  that  G(z)*  *  0;  thus  by  the  definition  of 

G(z)  in  (2.6),  ♦  must  first  satisfy  the  basic  scheme  associated  with  the 

resolvent  equation  g 

(2.30)  (Q_i  -  Tl 

x  o-0 


*  0,  v  *  1,2,3,. . . 


-  35 


Equation  (2.30)  is  an  ordinary  difference  equation  with  constant  coefficients; 
hence,  the  most  general  form  of  an  eigenvector  4>  satisfying  (2.30)  is  given  by 


(2.31) 


£=1  k=0 


£,k‘  £,k 


(v)< 


v  *  -r+1. 


Here,  =  k^(z),  are  the  distinct  inner  solutions  of  the  characteristic  equation 
(2.10)  each  with  corresponding  multiplicity  m  =  m  (z);.  P  (v)  are  arbitrary 

X  JC  X  ,K 

polynomials  in  v  with  deg[P  (v)]  =  k;  and  c  are  free  parameters  to  be 

1C  y  K  X  y  K 

determined,  where  by  Lemma  2.1,  their  precise  number  is 

n 

=  r- 

1=1 


REMARK  2.3'  We  note  that  the  splitting  property  mentioned  in  Lemma  2.1,  implies 
that  for  |z|  >  x,  the  inner  solutions,  <^(z),  satisfy  the  strict  inequality 
j ( z )  |  <  1,  hence,  the  eigenvector  $  given  in  (2.31)  is  in  £2(x).  Thus, 
the  existence  of  a  generalized  eigenvalue  z,  is  possible  only  for  z  lying  on 
the  unit  circle,  jz|  =1.  Furthermore,  operating  under  hypothesis  (2.11),  this 

possibility  is  reduced  in  the  dissipative  case,  to  the  single  point,  z  *  1,  since 

* 

the  splitting  property  in  this  case  is  maintained  also  for  |z|  =1,  z  f  1. 

We  now  make  a  particular  choice  of  the  polynomials  k(v)  in  (2.31), 
which  later  proved  useful.  We  choose 


* 


so  that  the  most  general  solution  of  (2.30)  which  is  used  as  eigenvector  of  the 
approximation,  is  of  the  form 


I 
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m.  -1 
n  1 


[2.32) 


-ZX!«.  *  »-  -r+l. 


k^o  i,k  k  1 


To  determine  the  parameters  a  ,  we  recall  that  being  an  eigenvector  of 
G,  <t>  must  also  satisfy  the  boundary  conditions 


(2.33) 


(G(z)$)  =  0  ,  y  =  0,-1,. . . , -r+l. 


The  operator  G(z)  is  defined  with  the  aid  of  the  operators  G^,  whose  operation 
at  the  boundary  points  is  given  by  (1.12);  hence  (2.33)  becomes 


( 2.  3*0 


T 

S  }u)  -  2Z  z'0-1  S^U))6  *  0,  u  -  0,-1 . -r+l. 

-1  o=0  a  v 


Inserting  (2.32)  in  (2,3*0,  we  finally  obtain 


m  -1 


(2.35) 


1=1  k=0  j=0 


0, 


u  »  0,-1,. . . ,-r+l, 


which  constitutes  a  linear  homogeneous  system  of  r  equations  in  the  r 
unknowns  o^  Clearly,  *  is  an  eigenvector  of  the  approximation,  if  and  only 
if  not  all  the  o^  ^  in  (2.35)  vanish,  that  is,  (2.35)  has  a  nontrivial 
solution. 

At  this  point,  we  associate  with  the  boundary  conditions  (1.12)  a  set  of 
rational  boundary-functions 


"fT> 


I 
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(2.36) 


R  (  Z,K  ) 

u 


=  ZJC 

J=0 


'  u ) 

J.-l 


T 


Z_z 

o=0 


Jo 


,  u  - 


0,-1 


,-r+l , 


which  are  uniquely  determined  by  the  boundary  coefficients  c 
the  scalar  functions 


Upon  defining 


(2.37a) 


c\v)kV*] 

3*0  Ja 


a  =  -1 ,0, . . .  ,t 


the  associated  boundary-functions  may  be  rewritten  as 


(2.3Tb) 


R  (z,<)  =  S(^(k)  -  5!z~0-1  S(u)(<),  u  =  0,-1 . -r+1. 

u  “  o=0  0 


Since 


j=o  1  a=o  J  '  K  ' 


the  system  (2.35)  is  cast  in  the  form 


m  -1 
n  l 


(2.38) 


E 


1=1  k=0  3k 


k  °t  k  =  °»  y  *  0,-1,. .. ,-r+l. 


It  follows  that  the  coefficient  matrix  of  this  system,  which  we  denote  by 


(2.39a) 


D  =  DU;^,...  vii^,...  ,mn) 


is  of  the  form 


I 
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D  =  [B(z,ic1,m1) ,  B(z,K2,m2),...,B(z,K  ,mn)], 
where  B(  z.k^  .tc^ ) ,  1  f  l  s  n  are  r*m£  dimensional  blocks  given  by 


(2.39 c) 


We  recall  that  *  is  an  eigenvector  of  the  approximation  if  and  only  if  (2.38) 
has  a  nontrivial  solution,  i.e.,  if  D  is  singular.  This  gives  us 
LEMMA  2.3. 

Let  z  with  |z|  z  1  be  given,  and  let  =  **.(*)»  1  S  l  S  n,  be  the 

corresponding  distinct  inner  solutions  of  the  characteri Stic  equation  (2.10) , 
each  with  multiplicity  m^  =  ( z ) . 

Then  z  is  an  eigenvalue  or  a  generalized  eigenvalue  of  the  approximation  if  and 
only  if 

det[D( z;x,  ,. . .  ,x  )]  ■  0. 

-L  n  i.  n 

By  Theorem  2.1,  the  stability  of  the  approximation  is  assured  if  and  only  if  it 
has  no  eigenvalues  nor  generalized  eigenvalues  z  with  |z|  *  1;  so  by  applying 


Lemma  2,3  we  finally  obtain  an  algebraic  formulation  of  the  stability  condition 
of  the  type  which  we  look  for. 

THEOREM  2.2. 

The  difference  approximation  (1.8) ,  (1.12)  o£  the  initial -boundary  value 
problem  ( 1 . T )  ij3_  3table,  if  and  only  if  for  every  z  | z |  i  1 ,  with  distinct  inner 
solutions  Kg ,  1  5  l  s  n,  each  with  multiplicity 

det[D(ziKlt. . . ,Kn;m1,. . . ,mn)]  4  0. 

Theorem  2.2  is  simplified  when  the  boundary  conditions  are  of  one-level  type, 
(see  (1.12b)),  i.e. , 

(2.40)  S^v^t+At)  =  0,  u  =  0,-1,. . .  ,-r+l. 

In  this  case  the  associated  boundary- functions  (2.37)  are 

r  (<)  =  s(;)(k)  = 

and  the  matrix  coefficient  D  in  (2.39)  is  given  by 


(2,41a, 


D  =  D(ic 


1  >  *  *  * 


*'cn’ml, 


,m  )  =  CB( ic  ,m ),...  ,B(k  ,m  )] 
n  n  n 


with 


I 


~  1*0  - 


(2. 1+loy 


B(VV 


RqU) 

VK) 

R_:(<) 

d 

■ 

’die 

L-r.l(‘L 

V"> 

d"'-1 

B.iU) 

*  v1 

die  1 

The  matrix  D  in  (2.1*1)  no  longer  depends  explicitly  on 
solutions  ic^  =  ic^z)  and  their  multiplicities  *  m^Cz), 
becomes 

COROLLARY  2.1. 

The  difference  approximation  (1.8),  ( 1 . 12b )  is_  stable 
every  z,  |z|  i  1  with  distinct  inner  solutions  1st: 

multiplicity  m^ 

det[D(ie^, . . .  .ic^;  m^,...,mn)]  4  0 


, ,  but  via  the  inner 
hence  Theorem  2.2 

if  and  only  if  for 
n,  each  with 


3.  SCHEME- INDEPENDENT  STABILITY  CRITERIA 


3.1.  Translatory  boundary  conditions  — ohe  determinant  condition 

In  the  previous  chapters,  we  dealt  with  boundary  conditions  of  the  general 
form  (1.12),  where  each  boundary  value,  v^(t+At),  is  determined  by  a  linear 
combination  of  computed  grid  values  which  is  dependent  on  the  position  of  the 
boundary  value  to  be  computed. 

In  this  chapter  we  start  discussing  translatory  boundary  conditions;  that  is, 
the  same  linear  combination  is  used  to  compute  the  boundary  values 
v^(t+At) ,  u  =  0,-1, ... ,-r+l,  independently  of  their  position.  In  other  words, 
the  translatory  case  is  characterized  by  applying  a  repeated  procedure,  where 
the  computation  is  done  by  translating  the  same  linear  combination  in  the  usual 
specified  order,  namely,  u  3  0,-1,. .. ,-r+l.  The  translatory  boundary  conditions 
are  thus  cast  in  the  form 

T 

S_1v  (t+At)  =  Z  S  v  (t-aAt),  y  ■  0,-1,. ..  ,-r+l,  t  i  TAt 
u  o=0  a  y 

(3'1)  V"  > 

. T 


where  the  coefficients  Cjo  are  no  longer  dependent  on  y. 

We  note  that  when  the  discrete  boundary  domain  is  reduced  to  the  single  point 
Xq,  such  as  in  the  case  of  the  widely  used  3-point  schemes  (i.e.,  r  *  p  *  1),  the 
computation  at  the  boundary  is  of  translatory  type  by  definition. 

Hereafter,  we  concentrate  on  searching  for  conditions  assuring  stability  in 
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the  translatory  case. 

The  rational  boundary- functions  associated  with  (3.1)  are 


(3.2a) 


q 

(z,k)  =  ~  c.(z)kv*^  tv  =  0,-1,. . . ,-r+l, 

w  0  J 


where  the  scalar  functions  Cj(z)  are  given  *>y 


(3.2b) 


JU>  "S'””1  *  °3,-l  ' 


o  j  J  a- 


In  particular,  for  p  =  0  we  get 


(3.3a) 


q.  t 

*n(z,tc)  *  c,(z)ie^  =  S  (ie)  -  ZZ  z~0_1S  (k)  , 
J=0  3  '’-0 


where  the  scalar  functions  ( k  ) ,  are  given  by 


q 


(3.3b) 


,«>-z 


.3 


'  jsj  V  ’  •  *  -1*0 . ’• 


By  (3-2a)  and  (3.3a),  we  have 


R^(z,ic)  =  KyR0(z,K)  ,  p  *  0,-1, . ...  ,-r+l. 


so  the  r  x  r  matrix 


D  —  D(z,k^,.  . .  ,<n , . . .  )  5  CB(z ) ,. .  .,B(t ) ) 


of  (2,39)  is  given  by  the  r*m  dimensional  blocks 


EU.Vm*) 


Rq(z.k) 

K  XR  (z,k) 


-r+l_  ,  % 

«  R0(z,<) 


RqU.O 
K  1Rq(z,k) 


K~r+1R0(z,te) 


Rq(z,ic) 
.  4 


’’  m-1 

3,  1 


k  I"*’1R()(z,ie) 


,lj4sn. 


The  fact  that  D  is  determined  now  by  the  single  boundary-function  Rq(z,k) 
enables  us  to  significantly  simplify  the  stability  condition  given  in  Theorem 
2,2  by  replacing  its  determinantal  criterion  with  the  following  scalar  condition. 
THEOREM  3.1. 

The  difference  approximation  (1.8),  (3.1)  of  the  initial -boundary  value 
problem  (1.7)  is.  stable,  if  and  only  if  for  every  z ,  | z |  >,  1,  with  corresponding 
distinct  inner  solutions  tc. ,  1  $  4  $  n,  we.  have 


(3.U) 


R0^  Z,<4 ^  ^  ^  *  l,2,.,.,n. 


PROOF.  Suppose  there  exists  zQ,  |zQ|  >.  1,  with  corresponding  inner  solution 

k  =  <.  (z),  j ic  |  *  1  violating  (3.^)  by  satisfying 
0  0  0 


Vz0»Kt  )  *  °* 

0 


Then  the  left  column  of  the  block  B{zn,ie.  ,mf  )  is  identically  zero,  hence  the 

0  *0  l0 


-  kh 


matrix  D  is  singular  and  by  Theorem  2.2,  the  approximation  is  unstable. 
Conversely,  suppose  (3.M  is  valid,  and  we  want  to  prove  stability, 
where  by  Theorem  2,2  it  suffices  to  show  that  for  every  a,  |z|  >,  1,  with 
distinct  inner  solutions  1  $  l  $  n,  each  with  multiplicity  m^,l<tfn, 

we  have 


For  that  purpose,  let 


det[D(z;x  ....  ,x  ;m  * 

i  n  i  n 


(3.5) 


u=-r+l 


®  »K-^) 

^-CxJr0(z,Ki)] 


»n-i 

X~^CK^Ro(z,Kn)3 

n 

_  n 


be  a  vanishing  linear  combination  of  the  rows  of  D.  The  vector  relation  in  (3.5) 
consists  of  the  r  scalar  equations 


V —  jk 

2_  <*  - k  CkJRoCz.ic^)]  -  OiO{k{  »t-l.  1  <  i  *  n. 


y=-r+l  3x 


which  we  write  as 


— Y  y  a  xjp11"1]  ■  Oj  0  (  k  <  m^-1,  1  i  t  $  n. 


u“-r+l 


(3.6) 
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By  our  hypothesis  (3.1*),  the  left  member  in  the  above  brackets  satisfies 


_  "I 

<l  nQ(z,Kt)  ¥  0  ,  1  *  t  $  n. 


Thus,  expanding  by  Leibniz'  rule  and  using  induction  on  k  >,  0,  we  find  that 
the  right  member  in  (3.6)  has  vanishing  derivatives  at  <  *  i.e.. 


=  0  ;  0  «  k  t  ».-l,  1  (  t  1  n. 

die  u=-r+l 


K~K. 


We  conclude  that  the  polynomial 


T(k)  *  ^  o  k 

U*-r+l 


r+p-1 


which  is  of  degree  r-1  at  most,  has  r  roots;  1  $  i  $  n  each  with 

multiplicity  m^.  Hence,  T(k)  i  0  and  the  coefficients  must  vanish.  By 

(3.5),  therefore,  the  rows  of  D  are  linearly  independent,  so  the  matrix  is 
nonsingular  and  stability  follows  by  Theorem  2.2. 

As  was  realized  in  the  previous  section,  if  (3.1)  is  reduced  to  the  one- 
level  case 


(3.7) 


S  ^v^(t+At) 


t  ca  , 

J*0  1 


j (t+At )  «  0,  y  «  0,-1,. . . ,-r+l. 


then  the  associated  boundary-functions 


R  00  *  c«  =  kWs  0*0,  U  *  0,-1,.. .  ,-r+l, 

u  _i 


r 
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cease  to  depend  explicitly  on  z,  and  Theorem  3.1  provide  us  with 
COROLLARY  3.1. 

The  difference  approximation  (1.8),  ( 3  •  7 )  ,is_  stable ,  if  and  only  if  for  every 
z,  I z |  i  1,  with  corresponding  distinct  inner  solutions  < 1  s  l  s  n,  we  have 

RqU*)  =  s_i(,c£)  *  0.  1  -  1,2,. ...n. 

Finally,  we  note  that  in  the  case  r  =  1,  which  was  mentioned  in  the  beginning 
of  this  section  as  translatory  one  by  definition,  the  matrix  representation  of  the 
results  which  obtained  in  Theorem  2.2  and  Corollary  2.2,  reduced  respectively 
to  the  scalar  results  given  in  Theorem  3.1  and  Corollary  3.1.  Indeed,  in  this 


case,  the  matrix  D  is  the  scalar  boundary- function  R^ . 


I 


-  UT  - 

3.2.  Scheme-independent  stability  criteria 

The  stability  criterion  given  in  Theorem  3.1,  involves  the  both  parts  which 
constitute  the  approximation  (1.8),  (3.1);  these  are  the  translatory  boundary 
conditions  (3.1)  which  generate  the  boundary-function  Rq(z,k),  and  the  (basic) 
scheme  (1.8)  which  induces  the  characteristic  equation  (2.10)  whose  z-dependent 
solutions  =  x  (z)  are  used  as  test  points  for  the  stability  of  the  approximation. 

Our  main  aim  in  this  section  is  to  provide  stability  criteria  which  do  not 
take  into  account  the  (basic)  scheme  (1.8),  but  instead,  are  given  solely  in  terms 
of  the  boundary  conditions .  In  such  a  way,  we  shall  be  able  to  answer  the  question 
whether  a  given  boundary  treatment  violates  the  stability  of  any  basic  scheme 
(as  an  approximation  to  the  pure-initial  problem-Assumption  III). 

REMARK  3.1.  We  emphasize  that  in  the  scheme-independent  stability  analysis 
carried  out  below,  it  is  always  assumed  that  the  (basic)  schemes  considered  obey 
the  four  basic  assumptions ,  Assumptions  I-IV,  which  were  originally  made  in  [6]. 

The  larger  part  of  this  section  discusses  the  (somewhat  simpler)  outflow  case. 

We  recall  that  these  are  the  outflow  scalar  components,  the  stability  of  which 
we  have  to  look  for  in  order  to  assure  the  stability  of  the  entire  vector 
approximation.  The  end  of  this  section  is  devoted  to  the  inflow  case,  where 
the  results  follow  easily,  merely  by  updating  the  results  previously  obtained  for 
the  outflow  problem. 

We  start  by  recalling  Lemma  7  in  [8]  which  discusses  the  behavior  of  the 
inner  solutions  ic^z)  in  'the  neighbourhood  of  the  point  z  *  1. 

LEMMA  3.1.  (Lemma  7  in  [8]) 

Consider  the  consistent  (basic)  scheme  (1.8)  as_  an  approximation  to  the  outflow. 


I 
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a  >  0,  (inflov,  a  >  0)  problem  (1.7).  Then,  its  associated  characteristic  equation 
(2.10)  has  exactly  one  outer  ( inner)  solution  ic  =  k(z)  which  satisfies 
<(zsl)  =  1. 

PROOF.  The  consistency  condition  ( 1.10a)  implies  that  z  =  1  is  a  solution  of 
the  eigenvalue  problem  (2.Ub) 


^  a  (z)elJf* 
j=-r  J 


*-i,e»o 


-  0, 


and  by  part  (ii) 
i.e. , 

(3.8) 


of  Assumption  III  (see  (2.5)),  z  *  1  must  be  a  simple  solution. 


t 

.Is5-] 


a^(z*l)  j  0 


]'  =  — 


dz 


3 


By  the  consistency  conditions  (1.10a),  (1.10b)  and  by  (3.8)  we  have 


(3.9a) 

(3.98) 


P(z’K,|z=x=l  Vz)|z=l  =  °* 

J=-r 

i  *  f  3*ju)U-i  *  *  °- 


j*-r 


Hence,  we  may  apply  the  Implicit  Function  Theorem  obtaining  that  in  the 
neighbourhood  of  z  =  1,  the  characteristic  equation  (2.10)  can  be  uniquely 
solved  for  k  as  a  differential  function  of  z;  that  is,  there  exists  a  single 
root-function  of  (2.10),  k  =  k(z),  which  satisfies 


(3,10a) 


k(z»1)  •  1. 


Applying  the  consistency  condition  (1.10b)  once  more  (see  (3.9b)  yields 


3P,  N 

rH  z,<)  i  . 

die  Z=K=1 


3P/ 


"Xa‘^Tr  aj(z)U=i  =  'Xa*3l{z’<)|z=<=i  • 


hence  the  root-function  k(z)  determined  above  satisfies 


o.iob)  «•(.) |„1--|r/|^|  -i/i*.  c 

1  1  z=ie=l 

Combining  (3.10a),  (3.10b)  implies  that  for  z  =  1+6, 
small,  we  have 


=  a[ _ 1 

dz  ' 


6  >  0  sufficiently 


ie(z)  =  l+(Xa)"16  +  Cf(52). 


Hence,  the  inequality 


(3.11)  |k(z}|  >  1  ( | < ( z ) |  <  1)  ,  a  >  0  (a  <  0), 

holds  in  the  right  real  neighbourhood  of  z  ■  1,  and  since  the  basic  scheme 
(1.8)  is  assumed  to  satisfy  the  von  Neumann  conditions  then  by  Lemma  2.11  this 
inequality  holds  for  all  z  with  |zj  >1.  Thus  k(z)  is  an  outer  (inner) 
solution  according  to  the  positive  (negative)  sign  of  the  coefficient  a 
which  completes  the  proof. 

In  the  course  of  our  discussion  about  scheme-independent  stability  criteria, 
we  introduce  two  additional  assumptions  complementing  the  first  four  already  made. 
We  will  show  that  the  new  assumptions  are  necessary  for  stability  and  provide 


I 
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scheme-independent  algebraic  tests  to  verify  their  validity. 

To  introduce  the  first  new  assumption,  let  the  scalar  functions  a^(z)  aJid 
Cj(z)  be  as  in  (l.ll)  and  (3.2b),  respectively. 


The  z-function  A(z),  given  by 


(3.12a) 


0  = 


|  i  aj(z)|  +  I  il  c.(z)|  , 
J=-r  J  J=0  J 


which  may  be  rewritten  in  the  form 


(3.12b) 


A(z)  =  |P{z,«=l)|  +  |rq(z,<=1) J  , 


satisfies 

(3.12c)  Mzse1^)  #  0  ,  0  <  ]u>|  f  w. 


REMARK  3.2.  We  note  that  Assumption  V  is  not  scheme- independent ,  since  the 
function  A(z)  in  (3.12)  depends  on  both-on  the  coefficients  a^(z)  determined 
by  tne  basic  scheme  (1.8)  and  on  the  coefficients  c^(z)  determined  by  the 
boundary  conditions  (3.1).  However,  it  should  be  pointed  out  that  although 
Assumption  V  depends  on  both  parts  of  the  approximation,  its  validity  can  be 
assured  by  considering  only  one  of  these  two  parts.  That  is,  (3.12c)  is  valid 
if  either  the  scheme -dependent  condition 
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(3.13) 


P(  z=e 


*•.»!)  .  £ 


J=-r 


aj(z=ei4>)  *  0  ,  0  <  |ip|  (  it. 


or  the  boundary-dependent  condition 


(3.1*0 


Rn(z=eltp,ic=l)  =  c  (z=ei<P)  ^0  ,  0  <  |ip|  $  it, 

J=0  3 


holds-  In  particular,  one  may  use  (3.1*0  as  a  scheme-independent  test  to  verify 
the  validity  of  Assumption  V. 

Verifying  the  validity  of  Assumption  V  becomes  much  simpler  when  the  basic 
scheme  or  the  boundary  conditions  are  either  two-leveled  or  three-leveled,  i.e., 
when  s  or  t  obtain  the  values  0  or  1.  This  is  the  content  of  the  next  lemma. 
LEMMA  3-2. 

Each  one  of  the  following  four  conditions.  is  sufficient  for  Assumption  V  to.  hold: 
(i)  -  The  basic  scheme  (1.8)  is  two-leveled,  i.e.  ,  s  *  0. 

(li)  -  The  boundary  conditions  (3.1)  are  two-leveled,  i.e.  .  t  »  0,  and  are 
accurate  of  order  (at  least)  zero. 

(iii)  -  The  basic  scheme  (1.8)  is.  three-leveled,  i.e. ,  s  *  1,  and  in  addition  we 
have 

(3.15a)  P(  z= 


*-l,<al)  =  / _  a  ( z=-l)  j  0. 

j«-r  i 


(iv)  -  The  boundary  conditions  (3.1)  are  three-leveled,  i.e.  ,  t  =  1,  are 
accurate  of  order  (at  least)  zero  and  in  addition  we  have 


Rrt(z«-l,ie«l)  5  e,(z«-l)  i  0. 

”  J»0  ” 


(3.15b) 


Thus,  roughly  speaking,  Assumption  V  is  automatically  fulfilled  in  the 
two-level  case,  and  has  to  be  verified  at  the  single  point  z  =  -1  in  the  three- 
level  one 

PROOF.  As  explained  in  Remark  3-2  above,  each  one  of  the  two  conditions  -  either 

(3.13)  or  (3.14)  is  sufficient  for  Assumption  V  to  hold;  thus  the  result  of 

our  lemma  follows  directly  from  Lemma  2.2. 

Indeed,  in  cases  (i)  and  (iii)  Lemma  2,2  implies  that  condition  (2.24) 

or  equivalently  (3.13)  holds  and  hence  Assumption  V  is  valid.  In  the  remaining 

cases,  (ii)  and  (iv),  the  accuracy  hypothesis  of  the  boundary  conditions  enables 

us  to  follow  the  proof  of  Lemma  2.2  replacing  the  functions  a.(z)  by  c  (z) 

J  0 

and  obtaining  that  condition  (3.l4)  holds.  Hence,  Assumption  V  is  valid  also 
in  these  cases. 

We  turn  now  to  discuss  scheme-independent  stability  criteria  and  let  us  start 
by  studying  schemes  of  dissipative  type.  The  important  point  in  the  stability 
analysis  of  such  schemes  is  the  fact  that  generalized  eigenvalues  z  with  |z1  3  1 
may  exist  only  at  the  single  point  z  ■  1  (see  Remark  2.3).  Recalling  also 
Lemma  3.1,  we  find  that  for  | z  j  3  1  all  the  corresponding  inner  solutions 
*t(z)  are  lying  inside  the  unit  disc,  i.e.,  | ( * ) |  <  1.  Indeed,  this  argument 
is  the  basis  for  our  next  theorem  discussing  scheme -independent  stability  criteria 
for  any  dissipative  (basic)  scheme. 

THEOREM  3.2. 

Consider  the  basic  scheme  (1.8)  of  dissipative  type  together  with  translatory 
boundary  conditions  as  an  approximation  to  the  outflow  problem  (1,7) . 

For  one-level  boundary  conditions,  t  ■' -1 ,  we  have 
(I)  -the  difference  approximation  (1.8),  (3.7)  is  stable  if  for  every  k  with 
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0  -  |*  <  ;i. 


13.26) 


n0U 


4. 

J  =  0  L 


¥  o. 


For  multi-level  boundary  conditions .  t  >  -1 ,  we  have 
■III)  -  the_  difference  approximation  (1.8),  (3.1)  is  stable  for  every  z  with  |  z|  >1 
and  every  <  with  0  <  | k j  <  1 , 


(3.17) 


Vz 


»<)  =  XZ 


_  c,(z)te^  ¥  0. 

j=0  J 


PROOF.  Take  an  arbitrary  z  with  \z\  i  1,  and  let  ic^z),  0  <  | «c z )  |  s  1, 
be  any  corresponding  inner  solution,  so  that 


(3.18) 


PU.k,  )  =  a^z)*^  *  0. 


J=-r 


In  order  to  assure  stability,  it  suffices,  according  to  Theorem  3.1,  to  show 
that 


(3.19a) 


Pq( 2 »K £ )  ¥  0. 


In  particular,  concerning  the  one-level  case  t  *  -1,-  the  boundary-function  Rc 
does  no  longer  explicitly  depend  on  z  and)  sufficient  condition  (3.19a) 
is  cast  in  the  form  (see  Corollary  3.1). 


(3.19b) 


i  0. 


"P 


For  the  case  where  the  inner  solution  k^(z)  is  inside  the  unit  disc,  i.e. 
0  <  | < ^ ( z ) |  <  1,  (3.19a)  and  (3.19b)  follow  respectively  by  hypothesis  (3.16) 
and  (3.17). 

Let  us  consider  then  the  case  where  the  inner  solution  ic^z)  is  lying  on 
the  unit  circle,  i.e.,  c^z)  ■  e1^,  0  i  |c|  f  ir. 

Our  assumption  of  the  dissipativity  of  the  basic  scheme  means  that  the  z- 
values  which  satisfy  (3.18)  with  inner  solution  of  the  form  ic^  =  e*^, 

0  <  |  s  rr,  obey  the  inequality 

>  3.20)  |  z(ic4=ei^)  |  <1,  0  <  |t|  S  v. 

Hence,  the  only  possibility  to  satisfy  (3-18)  by  an  inner  solution  of  the  form 
and  by  z  with  |z|  i  1,  is  the  possibility  of  =  e^j^^  =  1, 

where  by  (3.20),  continuity  implies  that  the  corresponding  z-value  satisfies 
(3.21)  I’W— “lE-o’l  4 

We  therefore  conclude  that  it  remains  to  verify  (3.19)  ia  the  case  where  ict=l 
and  the  corresponding  z-values  are  lying  on  the  unit  circle,  z  *  e^, 

0  (  |<p|  (  ir. 

Now,  since  the  dissipative  scheme  (1.8)  is  consistent  with  the  outflow 
problem  (1.7) •  then  according  to  Lemma  3.1,  ■  1  is  excluded  as  an  inner 

solution  corresponding  to  *  ■  1.  For  the  remaining  z-values,  s»ei<<>,0<|<p|<n, 
which  may  be  taken  into  consideration,  we  have  by  (3.18)  and  by  Assumption  V 
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(3.22) 


:z=e1U,,«  =l)j 


A(z*e^)  4  0  ,  0  <  |ipj  s  *  . 


That  is. 
The 
provides 
Thus  for 


(3.19)  is  valid  also  in  this  case,  and  stability  follows, 
first  part  of  Theorem  3.2  discussing  one-level  boundary  conditions, 
a  relatively  simple  stability  criterion,  as  it  depends  on  one  variable, 
example ,  m-order  extrapolated  boundary  values 


< . 


(3.23)  ( I-E)m+1v^  *  0  ,  u  -  0,-1,. . . ,-r+l 

can  be  easily  checked  as  satisfying  (3.16),  since  we  have 


R0(k)  s  (l-te)®*1  4  0  ,  0  <  | k |  <  1. 


In  particular,  in  the  case  of  a  two-level  basic  scheme,  where  the  validity 
of  Assumption  V  follows  by  Lemma  3.2,  we  obtain  the  following  well-known 
result  [3],  [73. 

COROLLARY  3.2.  (Theorem  5.2  in  [7]). 

The  two-level  dissipative  scheme  (1,8)  together  with  extrapolated  boundary 
values  as  given  in  (3.23)  constitute  a  stable  approximation  to  the  out  flow 
problem  (1.7). 

In  the  last  corollary  we  required  that  the  basic  scheme  will  be  two-leveled 
in  order  tc  assure  the  validity  of  Assumption  V.  When  vs  turn  to  extend  Corollary 
3.2  to  the  general  multi-level  case,  we  find  that  Assumption  V  is  indeed  a 
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necessary  one,  as  shown  by  the  following  example: 

Consider  the  three-level  5-point  dissipative  basic  scheme 

(3. PM  v  (t+At)  =  [l-r|{E-I )2I I-E-1)2]v  (t-At)  +  Xa(E-E-1)v  (t), 
v  lb  v  v 

v  =  1,2,...,  Xa  s  1-C ,  €  <  1 , 


€  2  —1  2 

which  follows  by  adding  the  dissipative  term -^-g(i -1 )  (I-E  )  v^(t-At),  to  the 
usual  Leap-Frog  scheme,  see  [ll.  Section  9].  For  both  schemes,  the  associated 
characteristic  equation  has  exactly  one  inner  solution  k  =«  ie(z),  satisfying 
tc(z=-l)  =  1  (see  [6,  Lemma  6.2]).  Now,  when  the  scheme  is  complemented  by 
m-order  extrapolation  of  the  boundary  values  Vg(t),v  ^(t),  the  approximation 
becomes  unstable,  as  follows  from  Corollary  3.1,  since  we  have 

R0C*{z)]|z=-l  *  (l-c,“2l  =  °* 

This  instability  is  explained  by  the  fact  that  the  approximation  fails  to  satisfy 
Assumption  V.  Indeed,  in  the  general  case  of  one-level  boundary  conditions 
which  are  at  least  zero-order  accurate,  such  as  the  m-order  extrapolation,  we  hcv 

by  (1.10a) 


0. 


Hence,  Assumption  V,  which  becomes 
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A  (  Z ) 


t  0 


z 


e1<P,  0  <  |ip|  s  u  , 


is  violated  by  our  scheme  (3.2U)  at  the  point  z  =  -1  since 


A(  z) 


z=-l 


=  0  . 


We  note  that  by  Lemma  3.2,  z  =  -1  is  indeed  the  only  possible  point  for  a 
three-level  scheme  to  violate  Assumption  V. 

It  was  already  mentioned,  that  the  stability  criterion  (3.16)  provided  in 
Theorem  3,21  which  discusses  one-level  boundary  conditions,  is  relatively  a 
simple  one,  since  it  involves  only  one  variable  ,  k.  Concerning  the  wide 
family  of  multi-level  boundary  conditions,  the  stability  criterion  provided  by 
Theorem  3-2II  in  (3.17),  is  more  complicated  since  it  involves  two  independent 
variables,  z  and  >c. 

This  motivates  us  to  look  for  z-independent  alternatives  to  Theorem  3.2II,  for 
both,  dissipative  and  particularly  non-dissipative  schemes,  complemented  by 
multi-level  boundary  conditions.  This  matter  will  occupy  the  remainder  of  our 
disscusion  about  the  outflow  problem. 

We  start  by  introducing  the  boundary-scheme  associated  with  the  boundary 
conditions  (3.1) 


T 

S_1vv(t+At)  «  Z.„w.  (t-oAt),  v  *  0,±1,±2,..., 


o*0 


0  V 


•t 


8,  .  .  -  -1.0 . t. 


(3-25) 


which  is  generated  by  extending  the  definition  of  the  boundary  values  in  (3.1), 
v^,  v  =  0,-1,. . . ,-r+l,  to  all  grid  points  ,  -«  <  v  <  °°. 

Since  the  boundary  conditions  (3.1)  were  assumed  to  be  multi-leveled,  i.e., 
t  is  non-negative,  it  follows  that  the  boundary-scheme  (3.25)  is  well-defined 
as  a  difference  scheme,  whose  values  are  computed  by  advancing  in  the  direction 
of  the  time-axis . 

As  was  realized  in  Section  2. 2, the  splitting  property  described  in  lemma 
2.1  was  the  key  in  investigating  the  solvable  (-Assumption  II),  scheme  (1.8), 
and  we  would  like  this  result  to  be  applicable  also  for  the  boundary-scheme 
(3.25).  For  that  reason,  we  require  the  solvability  of  the  boundary-scheme 
by  making  the  following  analogy  of  Assumption  II. 

ASSUMPTION  VI. 

The  boundary-scheme  (3*25)  is  solvable;  that  is,  there  exists  a  constant 
Kq  >  0,  such  that  for  every  ¥  €  l^[x) ,  there  is  a  unique  solution,  w,  w  £  i-2(x) 
for 

(3.26a)  S-lWv  3  »  v  *  -r+1  ,-r+2, . . . , 

with 

(3.26b)  iwl  l  5  Kq..  'v  5  ^  . 


REMARK  3.3.  To  assure  Assumption  VI,  one  may  use  Lemma  1.1  which  characterizes 
solvability.  Recalling  the  notation  there,  the  index  rQ  in  (1.16)  equals  zero 
in  the  case  of  the  right-sided  boundary-scheme  (3.26a),  since  by  (1.13)  ve  have 
cQ  ^  y  0.  Thus,  we  may  apply  the  second  part  of  Lemma  1.1,  replacing  the  scalar 
function  Q_|(<)  in  (1.17)  by  S_^(r)  given  in  (3.3b),  to  obtain  that 
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(3.27) 


;,(<)  e  5Z  c,  *  0  »  o  <  M  * 

x  i-n  J  »  x 


is  a  sufficient  condition  for  the  solvability  of  the  boundary-scheme  (3.25)* 

In  particular,  solvability  follows  for  explicit  boundary  conditions,  where 
S  ^(k)  =  constant. 

Before  continuing  we  want  to  associate  with  the  boundary  conditions  (3.1) 
two  concepts  -  the  von  Neumann  condition  and  dissipativity  -  concepts  which 
were  previously  associated  with  the  basic  scheme  (1.8).  We  will  find  it  quite 
attractive  to  express  our  forthcoming  stability  criteria  in  terms  of  these 
well-understood  and  easily  checkable  concepts. 

The  boundary-scheme  (3.25)  has  the  associated  characteristic  equation 


(3.28) 


V*»K)  -  21 

J-0  J 


by  means  of  which,  satisfying  the  von  Neumann  condition  and  dissipativity  make 
sense.  Upon  linking  these  properties  with  the  boundary- conditions  we  get: 

The  boundary  conditions  (3.1)  are  said  to  satisfy  the  von  Neumann  condition 


(3.29) 


y2,K=eU)  =  c  (z)elK  #  0  ,  |a|  >  1,  0  $  Id  f  * . 


and  are  said  to  be  of  dissipative  type  if 


RQ(z,K=ei5)  =  c.lzje1^  i  0  ,  |z|  i  1,  0  <  |c|  * 
J»0  J 


(3.30) 
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We  recall  that  Lemma  2.1  is  valid  for  any  solvable  scheme.  Applying  that 
lemma  for  the  solvable  boundary  scheme  (3.25),  yields  the  following  result. 

LEMMA  3-3 

For  the  boundary  conditions  (3.1)  which  satisfy  the  von  Neumann  condition, 
we  have  for  z,  |z|  j  1, 

(3.31)  R0(z><)  =  leW  t  0  ,  0  <  (  k  |  «  1. 

J-0  J 

PROOF  Take  an  arbitrary  z,  |z|  >1  and  consider  the  (polynomial)  characteristic 
equation  (3,28).  By  Lemma  2.11  all  of  its  non-zero  solutions  ic  *  ie(z) 
satisfying  |tc(z)|  >  1.  Thus  there  are  no  non-zero  solutions  of  (3.28)  in  the 
closed  unit  disc,  i.e,,  (3.31)  holds. 

REMARK  3 . U >  In  the  course  of  proving  Lemma  2.1  we  used  Assumption  I,  according 

tc  which  a  ( z)  4  0  ,  ( z|  a  1 ,  r  denoting  the  number  of  left  spatial  mesh 

points  that  the  basic  scheme  rests  on.  As  explained  in  Remark  2.1,  this  condition 

is  required  in  order  to  assure  that  multiplying  the  charact eristic  equation  by 
x* 

factor  ic  does  not  yield  additional  zero  solutions  <  *  0.  We  note  that  upon 
applying  Lemma  2.1  for  the  right-aidad  boundary  scheme  (3.25)  as  done  in  Lemma  3.3, 
we  are  free  from  requiring,  analogously  to  Assumption  I,  that  cQ(z)  4  0,  |z|*l, 

x* 

since  the  index  r  is  vanished  in  this  case,  i.e.,  k  =  1. 

By  continuity  arguments,  Lemma  3.3  implies  the  following  immediate  result. 
COROLLARY  3.3. 

For  the  boundary  conditions  (3.1)  which  satisfy  the  von  Neumann  condition, 
we  have  for  z,  |z|  a  1 
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5  1>c  I 


Rn!z’K  '  -  2Z  r:,(s)KJ  *  0  ,  0  <  |k[  <  1. 


J=0 


j 


L-iQl’Z*  Lemma  3.3,  the  solutions,  k  =  k  (z),  of  the  characteristic  equatic 

J  J 

'  ' .  ?8  )  satisfy  for  z,  |z|  >  1,  |<  (z)|  >  1. 

J 

Hence,  for  z,  |z|  j  1,  these  continuous  solutions  satisfy  |  *  ,(z)|  >  1,  ar.d 

J 

the  result  follows 


Combining  the  last  corollary  with  Theorem  3.2,  we  obtain  the  following  scheme- 
independent  stability  criterion  of  the  desired  type. 

THEOREM  3.3- 

The  basic  scheme  (1.8)  of  dissipative  type ,  together  with  the  boundary 
conditions  (3.1)  which  satisfy  the  von  Neumann  condition,  constitute  a  stable 
approximation  to  the  outflow  problem  ( 1 . T ) . 

PROOF ,  Since  the  boundary  conditions  (3.1)  satisfy  the  von  Neumann  condition, 
then  by  Corollary  3.3  we  have  (3.32)  which  by  the  second  part  of  Theorem  3.2 
is  sufficient  for  stability. 

The  last  theorem  provides  a  scheme-independent  stability  criterion  for 
difference  approximations  whose  basic  scheme  is  limited  to  be  of  dissipative  type. 
We  turn  now  to  the  general  case,  of  basic  schemes  which  are  not  necessarily 
dissipative.  In  particular,  we  refer  to  the  case  of  non-dissipative  schemes 
where  unlike  the  dissipative  case  all  z  lying  on  the  unit  circle  may  serve  as 
generhlized  eigenvalues. 

THEOREM  3.  **. 

The  basic  scheme  (1.8)  together  ■■/ith  the  boundary  conditions  (3.1)  of 

dissipative  type,  constitute  a  stable  approximation  to  the  out  flow  problem  (1.7,, 


r 


Tr 


I 


I 


-  62  - 


PROOF  Take  an  arbitrary  z  with  |z|  >  l  and  let  <  [z) ,  0  <  |<^(z)|  f  1,  be 
any  corresponding  inner  solution,  so  that 


1 ?. 33) 


P(  z.tc^) 


t 

J— r 


0. 


In  order  to  assure  stability,  it  suffices,  according  to  Theorem  3.1,  to 
show  that 

^  q. 

(3-3M  Rr 

J=0 


*  °- 


We  first  note,  that  since  the  boundary  conditions  were  assumed  to  be  of 
dissipative  type,  they  particularly  satisfy  the  von  Neumann  condition;  so  Lemma 
3  3  and  Corollary  3.3  may  apply  to  our  case. 

Now,  for  z-values  outside  the  unit  disc,  |z|  >1,  (3.3^)  follows  from 

Lemma  3  3,  and  for  an  inner  solution  which  is  inside  the  unit  circle, 

0  <  1*^1  <  1,  (3-3M  follows  from  Corollary  3-3. 

Therefore,  it  remains  to  verify  (3.31*)  for  the  case  that  both  z  and  < 
are  lying  on  the  unit  circle,  i.e., 

(3.35)  z  =  eltP  ,  0  s  |<p|  $  it  ;  =  e1?  ,  0  j  |  £ |  s  * . 


For  an  inner  solution  of  the  form  k  =  e1^,  £  ¥  0,  (3.3^0  follows  from  the 
dissipativity  of  the  boundary  conditions  (see  (3.30)).  So  let  us  consider  an 
inner  solution  of  the  form  < ^  *  e^j^g  a  1. 

By  Lemma  3.1,  upon  approximating  the  outflow  problem  (1.7),  *  1  is 


vr 


1 
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excluded  as  an  inner  solution  corresponding  to  2=1;  and  for  the  remaining 
z-  values,  z  =  ,  0  <•  J  op]  5  r ,  which  may  be  taken  into  consideration,  we  have 

by  (3-33)  and  by  Assumption  V 

(3.36)  |R0(z=ellVjl=l)|  =  A(z=e143)  ?  0  ,  0  <  |cp|  n  . 

That  is,  (3.3M  is  valid  also  in  the  remaining  case,  and  stability  follows. 

Combining  Theorem  3-3  and  3.^  we  immediately  obtain  the  following  summary 

result . 

COROLLARY  3.^. 

Consider  the  basic  scheme  (1.8)  together  with  the  boundary  conditions  (3.1) 
which  satisfy  the  von  Neumann  condition ,  as  an  approximation  to  the  outflow 
problem  (17).  I_f  either  the  basic  scheme  (1.8)  or  the  boundary  scheme  (3.25)  is 
dissipative ,  then  the  approximation  is  stable . 

We  note  that  the  stability  properties  of  the  boundary  scheme,  namely, 
dissipativity  and  the  von  Neumann  condition,  are  often  known  in  advance .  Thus, 
in  applying  the  last  scheme-independent  stability  criteria  summarized  in  Corollary 
3.i*,  then  beside  the  four  basic  assumptions  (Assumptions  I-IV)  which  the 
approximation  is  always  assumed  to  satisfy  (see  Remark  3.1),  it  remains  to  verify 
the  validity  of  the  additional  assumptions,  Assumptions  V  and  VI. 

For  the  purpose  of  assuring  these  assumptions,  one  may  use  Lemma  3.2  and 
Remark  3.3  which  imply  in  particular  that  Assumption  V  is  automatically  fulfilled 
in  the  case  of  two- level  boundary  conditions,  and  that  the  solvability  assumption 
VI  i3  automatically  fulfilled  in  the  case  of  explicit  ones. 


J 
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Yet,  referring  to  the  general  multi-level  implicit  case,  then  Assumptions 
V  and  VI  are  indeed  necessary  for  the  validity  of  the  scheme-independent 
stability  criteria  given  in  Theorems  3.3  and  3.1*. 

Concerning  the  first  of  these  two,  then  by  Lemma  3.2  it  follows  that 
Assumption  V  is  automatically  fulfilled  in  the  two-level  case.  When  we  raise 
the  question  whether  stability  is  maintained  also  in  the  general  multi-level 
case  involving  more  than  two  time  steps,  the  following  example  shows  that  the 
answer  to  that  question  is  negative.  That  is,  Assumption  V  is  indeed  a 
necessary  one. 

Consider  the  non-dissipative  Leap-Frog  scheme 
(3.37a)  vy( t+At )  =  v^(t-At )  +  Xa* Cv^^Ct )-v^_1(t ) ]  ,  v  *  1,2,..., 

together  with  the  solvable  consistent  boundary  condition 
(3.3Tb)  vQ(t+At)  =  vQ(t-At)  +  2Aa- [^(t-At )-vQ(t-At ) 3 . 

The  boundary  function  associated  with  (3.37b)  is  given  by 
(3.38)  Rq( z ,k)  =  1  -  z_2*[l+2AaU-l)]  , 

and  its  z-roots,  z  =  z(k),  satisfy  for  0  <  Xa  5  0.5 


(3.39) 


|z2U*eU)!2<  (|l-2Aa|  ♦  |2Aa|)2  ■  1  ,  0  <  |£|  *  « 
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Thus,  the  explicit  (and  hence  solvable)  boundary  condition  (3.37b)  is  of 
dissipative  type;  so  approximation  (3.37)  fulfills  the  requirements  of  both 
Theorem  3.1*  and  Assumption  VI. 

At  the  same  time,  the  approximation  (3.37)  is  unstable.  Indeed  the 
characteristic  equation  associated  with  (3.37a) 

P(z,tc)  =>  l-z^-z'^XaU-*"1)  =  0, 

has  exactly  one  inner  root- funct ion ,  tc  «  x(z),  satisfying  x(z=-l)  =  1 
(see  [6,  Lemma  6.2]),  and  by  inserting  it  into  the  boundary-function  (3.38), 
we  get 

(3.1*0)  Rn(z=-i,K=l)  =  l-z"2-[l+2Xa(*:-l)],  =  0. 

0  Z=-l 

'<=1 

Thus,  the  approximation  (3.37),  is  unstable  since  it  violates  the  necessary 
stability  condition  (3-1*)  at  the  point  z  =  -1. 

The  instability  of  approximation  (3.37),  despite  that  it  fulfills 
the  requirements  of  both  Theorem  3.1*  and  Assumption  VI,  is  explained  by  its 
failure  to  satisfy  also  Assumption  V  as  follows  from  (3.38): 

(3.Ul)  a(z=eitP|q>:_n)  *  |Rq( z»-l,*c=l)  |  =  0. 

We  remark  that  according  to  Lemma  3.2,  the  only  possibility  of  the  three- level 
approximation  (3.37)  to  violate  Assumption  V  is  at  the  single  point  i  »  -1, 
as  we  have  indeed  found  in  (3.1*1). 


Concerning  the  solvability  Assumption  VI,  its  necessity  can  be  shown  by 
considering  any  two-level  3-point  dissipative  basic  scheme  together  with  zero 
order  accurate  boundary  condition  of  the  form 

(3.42)  Vq( t+At )  -  Bv^t+At)  ■  vQ(t)  -  Bv  (t)  ,  B  >  1. 

The  boundary- function  associated  with  (3.42)  which  is  given  by 
(3.1*3)  R0(z,*O  =  ( l-z-1)  •  ( 1-Bic)  , 

satisfies  RQ(  z.iese1^)  4  0  ,  | z |  >  1,  0  s  | £ |  $  *,  hence  by  (3.29)  the  boundary 
condition  (3.42)  satisfies  the  von  Neumann  condition  and  the  entire  approximation 
fulfills  the  requirements  of  Theorem  3-3.  Furthermore,  since  the  basic  scheme 
was  assumed  to  be  two-leveled,  then  by  Lemma  3.2,  Assumption  V  is  fulfilled  as  well. 
Yet,  the  approximation  is  unstable  since  the  boundary-function  (3.1*3)  vanishes  at 
z  =  1  independently  of  te-values;  hence  the  necessary  stability  condition  (3.1*) 
is  violated.  This  instability  is  explained  by  the  failure  of  the  boundary  scheme 
associated  with  (3.42)  to  be  solvable.  Indeed,  recalling  the  solvability  definition 
in  (3.26),  then  by  taking  V  =  0  in  (3.26a)  we  find  that  the  grid  function 

w  =  |b  VWq|v_  r+1£J.0(x)  with  arbitrary  Wq,  satisfies  S  ^wy  =  wv~8wv+1  =•  0,  vj-r+1. 

Thus,  we  have  neither  the  uniqueness  nor  the  boundedness  which  is  required  in 
(3.26b). 

Our  study  of  the  outflow  problem  is  completed,  and  we  turn  now  to  discuss 


I 


-  67  - 


some  remarks  concerning  the  inflow  one. 

We  first  note  that  all  the  results  which  were  discussed  in  previous  sections 
go  over  unchanged  except  for  Lemma  3.1.  The  result  of  this  lemma  discussing  the 
behavior  of  inner  solutions  in  the  neighbourhood  of  z  =  1,  depends  on  whether  we 
approximate  the  outflow  problem  or  the  inflow  one. 

In  the  outflow  case.  Lemma  3.1  is  used  to  exclude  the  possibility  of  ^  ■=  1 
to  serve  as  an  inner  solution  corresponding  to  z  *  1.  In  the  inflow  case, 
however,  the  situation  is  Just  the  contrary;  that  is,  according  to  Lemma  3.1,  z  =  1 
has  always  exactly  one  corresponding  inner  solution  k.  =  1. 

Now,  we  recall  that  all  our  stability,  criteria,  particularly  the  scheme- 
independent  ones,  were  obtained  by  applying  Theorem  3.1  which  characterizes 
stability  by  requiring  that  for  every  t  with  |z|  >,  1  and.  every  corresponding 
inner  solution  =  <t(z),.l  si*  n,  we  have  R0(z,k^)  #  0. 

We  therefore  conclude  that  when  dealing  with  the  inflow  problem,  all  our 
previous  stability,  criteria  still  hold  upon  making  the  additional  requirement 

R0(z»l,iei=l)  i  0, 

a  requirement  which  was  automatically  excluded  by  Lemma  3.1,  in  the  outflow  case. 

Thus  for  example,  referring  to  the  summary  result  in  Corollary  3.*+,  we  obtain 
for  thq  inflow  problem 
COROLLARY  3.5. 

Consider  the  basic  scheme  (1.8)  together  with  the  boundary  conditions  (3.1) 
which  satisfy  the  von  Neumann  condition,  as  an  approximation  to  the  inflow  problem 
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( 1 . T) .  If.  either  the  basic  scheme  (1.8)  or  the  boundary  scheme  ( 3.25 )  is 
dissipative  and  if  in  addition  ve  have 


(3.^2) 


Rq(z=1,k=1)  4  0, 


then  the  approximation  is  stable . 

We  note  that  when  the  additional  condition  (3. 1*2) 


R0(z»* ) 


Z=K=1 


cj(z) 


Z=1 


4  0 


holds,  then  according  to  ( 1.10a)  the  boundary  conditions  must  be  inconsistent 
and  in  fact  have  no  accuracy  with  respect  to  the  differential  equation.  This 
indeed  makes  sense  since  one  cannot  expect  the  stable  approximation  (1.8),  (3.1) 
whose  values  are  uniquely  determined  in  the  quarter-plane  x,  t  >,  0,  to  be 
consistent  with  the  inflow  problem  (1.7)  which  is  not  uniquely  determined  unless 
extra  boundary  data  is  supplied  as  given  in  (1.7b).  Thus  in  general,  consistent 
boundary  conditions  of  translatory  type  approximating  the  outflow  problem  are 
of  no  value  when  dealing  with  the  inflow  one.  Yet,  there  is  one  important  case 
which  we  shall  now  consider.  That  is,  when  the  missing  boundary  values  are 
computed  via  summations  of  the  form  (see  [7,  Theorem  6]) 


(3.U3) 


gu(t)et2(t),  W  =  0,-1,. . . ,-r+l. 


For  the  purpose  of  determining  stability,  we  consider  the  boundary- function 
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q. 

associated  with  the  homogeneous  part  of  (3.1+3),  which  is  given  by  R^k)  = 

We  have  R0(>c)  ^  0  for  all  tc  ^  0  and  hence  by  Corollary  3.1,  the  stability 
of  both  the  outflow  and  particularly  the  inflow  approximation  (1.8),  (3.^3) 
is  assured.  Setting  q  to  be  zero,  we  obtain  the  well-known  result  C T ] ,  C 10 3 
of  the  unconditional  stability  of  the  approximation  whose  boundary  values  are 
determined  by  arbitrary  bounded  inhomogeneous  terms  v^(t)  a  gy(t)  * 
u  =  0,-1,. . . ,-r+l.  Indeed,  this  result  was  mentioned  earlier  in  Section  1.1, 
where  it  was  used  to  assure  the  unconditional  stability  of  the  inflow  scalar 
components  of  the  vector  approximation  (1.5),  (1.3). 
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4.  EXAMPLES  OF  SCHEME- INDEPENDENT  STABILITY  INVESTIGATIONS. 

In  this  chapter  we  study  some  examples  of  translatory  boundary  conditions 
which  together  with  corresponding  basic  schemes  constitute  stable  approximations 
to  the  outflow  problem  (1.7). 

For  that  purpose,  we  apply  the  scheme-independent  criteria  of  the  previous 
chapter,  so  that  stability  is  not  restricted  to  a  specific  basic  scheme.  That 
is,  the  acquired  stability  is  valid  for  a  family  of  approximations  which  consists 
of  the  boundary  conditions  together  with  any  basic  scheme  having  some  general 
property  (the  "familial"  property)  such  as  dissipativity ,  two-levelness,  etc. 
REMARK  4.1.  It  is  of  course  understood  that  beside  requiring  the  basic  schemes 
to  satisfy  some  general  ("familial")  property  which  follows  from  the  scheme- 
independent  stability  analysis,  all  basic  schemes  considered  must  satisfy  the 
four  basic  assumptions.  Assumptions  I-IV  (see  Remark  3.1).  In  particular,  we 
refer  to  the  stability  assumption,  Assumption  III,  which  may  lead  to  impose  some 
restriction  on  the  time  step  used  At. 

We  note  that  the  stability  criteria  given  in  Theorems  3.3  and  3.4  are 
independent  of  the  index  r,  which  denotes  the  number  of  boundary  values  to  be 
computed  at  each  time  level  v^(t),  y  =  0,-1,. . . ,-r+l.  Hence,  verifying  the 
stability  in  the  simpler  case  of  computing  a  single  boundary  value,  v^(t), 
(complementing  for  example,  a  3-point  basic  scheme)  requires  no  more  effort  than 
the  stability  verification  in  the  general  translatory  case  of  r  boundary  values, 
>  1,  complementing  a  basic  scheme  of  the  general  form  (1.8). 

The  above  observations  are  particularly  relevant  to  those  boundary  treatments 


which  are  considered  below,  and  whose  stability  is  already  discussed  in  the 


literature.  However,  the  verification  of  stability  given  below  has  two  specific 
features:  first,  because  the  stability  investigation  is  independent  of  the  basic 
scheme  and  of  the  solutions  of  the  corresponding  characteristic  equation,  then 
the  procedure  becomes  much  shorter;  secondly,  the  translatory  nature  of  the 
boundary  treatment  assures  that  the  acquired  stability  is  valid  for  any  basic 
stable  scheme  and  not  necessarily  for  3-point  ones. 

Let  us  turn  then  to  the  examples  themselves,  and  consider  first  boundary 
conditions  which  complement  any  dissipative  basic  scheme. 

EXAMPLE  4.1.  (example  (6.11)  in  [6]).  Let  the  boundary  conditions  be  determined 
by  oblique  Lagrangian  extrapolation  of  order  m-1: 


(4.1a) 


v^(t+At) 


XU 

-z 

i=l 


(“)(~l)J+1Vu+J[t-(J-l)At] 


u  =  0,-1,. 


-r+1. 


The  boundary-function  associated  with  (4.1a)  is  given  by 


(4.1b)  R  (z,x)  =  1-  1)  (®)(-l)J+1z"’3icJ  =  (l-z_1ic)m, 

J«1  J 


and  by  equating  to  zero,  we  get  that  the  z-solutions  of  (4.1b)  satisfy 


|zU*eU)|  »  |eU|  =  1  ,  Vt. 


Thus,  the  explicit  and  hence  (by  Remark  3.3)  solvable  boundary  conditions  (4.1a) 
are  of  non-dissipative  type  so  they  satisfy  the  von  Neumann  condition.  According 
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to  (4.1b)  we  also  have  z  ( <=1 )  =  1,  so  A(z=e1(^)  z  |  R^(  z=e1^>,<=l )  |  >0 ,  0<|<p|sir, 
i.e..  Assumption  V  is  valid  and  stability  follows  by  Theorem  3.3. 

EXAMPLE  4.2.  (example  (6.3c)  in  [6],  example  (3.4)  in  [l4]).  Let  the  boundary 
conditions  be  generated  by  the  Box-Scheme,  i.e.. 


(4.2a)  v^( t+At )  +  v^+1(t+At)  -  Xa*[v^+^(t+At)  -  v^(t+At)]  = 

=  v^(t)  +  vp+1('t)  +  Xa*[v^+1(t)-v^(t)]  ,  u  =  0,-1,. .. ,-r+l. 
The  boundary- function  associated  with  (4.2a)  is  given  by 
(4.2b)  Rq(z,k)  =  l+K-Xa(x-l)  -  z_1*[l+ic+Xa(ic-l)]  , 


-  by  equating  to  zero,  we  get  that  its  z-solutions  satisfy 


z(ic=e1^ )  | 


l+elg+Aa(eU-l) 

l+ei5-Xa(ei5-l) 


1  ,  ve. 


Thus,  boundary  conditions  (4.2a)  are  of  non-dissipative  type  so  they  satisfy 
the  von  Neumann  condition.  Approximating  the  outflow  problem  (a>0)  we  have 


Re[S_1(K)]  ■  1+Re(<)  ♦  Xa*[l-Re(ie)]  ¥  0  ,  |ie|  f  1. 

Hence  (see  Remark  3.3),  the  boundary  conditions  (4.2a)  are  solvable  (i.e.. 


Assumption  VI  holds),  and  since  they  are  also  two-leveled  then  (by  Lemma  3.2) 


I 
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Assumption  V  holds  as  well.  Therefore  we  may  apply  Theorem  3.3  obtaining 
stability. 

In  view  of  the  stability  discussion  in  Examples  4.1,  4.2,  it  follows  that 
any  dissipative  basic  scheme  together  with  translator y  boundary  conditions  which 
are  generated  by  oblique  extrapolation  or  by  the  Box-Scheme ,  constitute  a  stable 
approximation  to  the  outflow  problem  (1.7). 

EXAMPLE  4 . 3 .  ([2]  ,  example  (6.2b)  in  [ 15 3 ) •  Let  the  boundary  conditions  be 
generated  by  the  right-sided  wieghted  Euler  scheme,  i.e., 

v  (t+At)  =  v  (t-At)  +  2Xa*[v  . (t  )-0. 5*  ( v  (t+At)  +  v  (t-At))], 
u  u  u+1  p  u 

U  —  0,-1,. . . , — r+ 1 . 


The  boundary-function  associated  with  (4.3a)  is 


(4.3b)  Rq( z  ,k  )  =*  l-z_2-2Xa*[K*z_1-0.5*(l+z  2)3» 


and  by  equating  to  zero,  we  find  that  its  z-solutions  are  given  by 


(4,3c)  zU'e1^)  =  ei?  ,  b(0  =  /(Xa)2+e-2i^[l-(  Xa)2]. 


In  order  to  assure  stability,  we  restrict  the  time-step  At  by  requiring 


the  Courant-Friedrichs-Levi  (CFL)  condition 
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(4, 3d)  0  <  Xa  s  1. 

We  note  that  the  CFL  requirement  (4.3d)  is  automatically  fulfilled  in  the  case 
of  an  explicit  basic  scheme,  since  by  our  Assumption  III,  the  basic  scheme  must 
satisfy  the  von  Neumann  condition. 

Having  (4.3d),  it  follows  that  |b(£)|  S  d»  0  $  |?|  $  w.  Thus,  the  explicit 
and  hence  (by  Remark  3.3)  solvable  boundary  conditions  (4.3a)  satisfy  the  von 
Neumann  condition 

|  z(<=elC)  I.  *  HKOl  <  !  ,  yg. 

According  to  (4.3c)  we  have  z(x=l)  *  4  eiip,  0  <  jip|  5  *,  so 

A(z=eiiP)  2  | RQ( z^e^.ie^l)  |  >  0,  0  <  |<pj  5  tv,  i.e.,  Assumption  V  is  valid.  We 
may  apply  now  Theorem  3.3  to  conclude  that  any  approximation  to  the  outflow 
problem  (1.7)  which  satisfies  the  CFL  condition  (4.3d)^\  consisting  of  a 
dissipative  basic  scheme  together  with  the  translatory  boundary  conditions  (4.3a), 
is  stable. 

We  turn  now  to  study  general  difference  approximation  consisting  of  any  basic 
scheme  (which  is  not  necessarily  of  dissipative  type)  together  with  the  following 
boundary  treatments. 

EXAMPLE  4.4.  (example  (6.3b)  in  [6]  ,  example  (3.2)  in  [l4]).  Let  the  boundary 
conditions  be  generated  by  the  right-3ided  explicit  Euler  scheme,  i.e., 

(4.4a)  vy(t+At)  =  v^(t)  +  Xa*[vy+1(t)-v  (t)3,  y  *  0,-1, . . .  ,-r+l. 

Tl) A  further  restriction  on  the  time  step  At  may  arise  from  the  stability 
requirement  made  in  Assumption  III  (see  Remark  4.1). 
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The  boundary-function  associated  with  (4.4a)  is  given  by 
(  U .  Ub )  Rq(  z ,k  )  =  l-2_1*[l+Xa(ic:-l)]  , 

and  by  equating  to  zero,  we  get  that  its  z-solutions  satisfy 

|z(<=e1?) |2  =  (1-Xa+Xa-cosC)2 +  (Xa-sin^)2 <  ( |l-Xa|+|Xa| )2  ,  0<  |c|  (*■ 
Thus,  requiring  the  CFL  condition 

(4.4c)  0  <  Xa  $  1, 

we  see  that  the  boundary  conditions  (4.4a)  are  of  dissipative  type.  Since  the 
boundary  conditions  (4.4a)  are  also  explicit  and  two-leveled  then  by  Lemma  3-2 
and  Remark  3-3,  both  Assumptions  V  and  VI  are  valid  and  stability  follows  £ppm 
Theorem  3.4. 

EXAMPLE  4 , $ .  (example  (3.3)  in  [l4]).  Let  the  boundary  conditions  be  generated 
by  the  right-sided  implicit  Euler  scheme,  i.e., 

(4.5a)  v^(t+At)  -  Xa*[v^+1(t+At )  -  v^(t+At)]  =  vy(t)  ,  u  =  0,-1,. . . ,-r+l. 

The  boundary-function  associated  with  (4.5a)  is  given  by 

Rq(  z  ,ic  )  «  l-XaU-D-z”1, 


(U.5b) 
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and  by  equating  to  zero,  we  get  that  its  z-solutions  satisfy  in  the  outflow  case 
(a>0) 

|  z(K=e'*'^ )  |  ^  =  [  ( 1+Xa-Xa  •  cos£) "  +  ( Xu*  sin£  )^3  ^  <  [|l+Xa| -|  Xa|  D  =  1,  0  <  |s|  f  it. 

Thus,  boundary  conditions  (*i.5a)  are  of  dissipative  type.  Approximating  the 
outflow  problem  (a>0)  we  have 

Re[S_1  ( k ) ]  =  1+Xa- l-Re(ic)]  4  0  ,  jx|  <:  1. 

Hence  (see  Remark  3-3),  the  boundary  conditions  (L.5a)  are  solvable  (i.e.. 
Assumption  VI  holds),  and  since  they  are  also  two-leveled  then  (by  Lemma  3.2) 
Assumption  V  holds  as  well.  Therefore  we  may  apply  Theorem  3.1*  obtaining 
stability. 

In  view  of  the  stability  discussion  in  Examples  h.L  and  U.  5,  we  may  conclude 
that  if  the  boundary  conditions  (3.1)  are  generated  by  a  stable  right-sided 
explicit  Euler  scheme  or  by  the  right-sided  implicit  Euler  scheme,  then  the 
entire  approximation  (1.8),  (3.1)  is  stable. 

EXAMPLE  h,6.  Let  the  boundary  conditions  be  of  the  form 

(U.6a)  (l+Xa)v^(t+At)  +  (l-Xa)v^+^(t+At)  = 

=  2Xa*[v  ..(t)-v  (t)3  +  (l-Xa)v  (t-At)  +  (l+Xa)v  ..(t-At)*  y  *  0,-1,. . . ,-r+l. 
y+1  y  y  y+1 


I 
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The  boundary- function  associated  with  (4.6a)  is 

(4.6b)  Rq(z,k)  =  1+Xa+ ( 1-Xa )k-z  ^  .2Xa(ic-l)-z  ^'[(1-Xa)  +  (l+Xa)*:], 


and  by  equating  to  zero,  we  find  that  its  z-solutions  are  given  by 


(4.6c) 


/  \ _  Xal • sin^/2  ±  cosg/2 

1  -Aai-ainE/2  +  cosC/2  ’ 


and  hence  (z(»c=ei^)|  =  1.  Thus,  the  boundary  conditions  (4.6a)  are  of 
non-dissipative  type  so  they  satisfy  the  von  Neumann  condition.  Furthermore, 
we  have  in  the  outflow  case  (a>0) 


Re[S  (*)]  =  1+Re(i«)  ♦  Xa*[l-Re(  ic )  ]  ^  0  ,  |k|  si, 

and  hence  (see  Remark  3.3)  boundary  conditions  are  also  solvable,  i.e..  Assumption 
vi  holds.  To  assure  stability  via  Theorem  3.3.  it  then  remains  to  verify 

Assumption  V. 

Now,  since  by  (4.6c)  we  have  z(k=1)  »  ±1,  hence 

A(z«eilp)  i  jR0U*>ei<P,K-l)  |  >  0  ,  0  <  |<p|  <  n  , 
and  therefore  the  validity  of  Assumption  V  follows  upon  requiring 


A(z=-1)  ■  |p(z=-l,te»l)|  i  0. 


I 
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Here,  P(z,k)  de:.'-tcs  as  usual,  the  characteristic  function  associated  with  the 
basic  scheme. 

Applying  Theorem  3.3  we  conclude  that  the  boundary  conditions  (4.6a),  which 
complement  any  dissipative  basic  scheme  whose  characteristic  function  P(z,k) 
satis fi es 

(4.6d)  P(z=-1,k=1) 


constitute  a  stable  approximation  to  the  outflow  problem  (1.7). 

Since  by  Lemma  2.2J  condition  (4.6d)  which  may  be  rewritten  in  the  form 


t 

j*-r 


aj(z=-l)  i  o 


is  automatically  fulfilled  in  the  two-level  case,  we  obtain  that  stability 
follows  whenever  the  boundary  conditions  (4.6a)  complement  any  two-level 
dissipative  basic  scheme.  We  note,  however,  that  in  the  general  case  of  multi¬ 
level  basic  schemes  involving  more  than  two  time-steps,  the  additional  requirement 
(4.6d)  is  indeed  a  necessary  one  as  shown  by  the  following  example: 

Consider  the  three-leveled  dissipative  scheme  (see  (3.24)) 


(4.7a) 


\ v ( t+At )  =  [I  -  y|^E-I )2(I-E~1)2]vv(t-At)  +  Xa(E-E_1)vv(t) , 


v  =  1,2,...,  Xa  s  1-C,  £  <  1, 


1 


and  let  the  boundary  value  Vg(t+At),v  ^(t+At)  be  computed  by  (4.6a),  i.e., 


I 
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(i*.7 


{1+Aa)v  (t+At) 

p 


(l-Aa)v  ( t+At '  = 
P+1 


=  2Aa' 


Cv  .. (t)-v  (t)]  + 

P+1  P 


[l-Aa)v  (t-At)  ♦  (1+Aa)v  ,(t-At),  p  =  0,-1. 

p  p+1 


The  characteristic  equation  associated  with  the  basic  scheme  (U.7a),  has 
exactly  one  inner  root  function  <  =  x(z)  satisfying  x(z=-l)  =  1,  and  by 
inserting  into  the  associated  boundary-function  we  get 


R0(j~1,k=1)  =  l+Aa-(l-Aa)x-z"1*2Aa(1c-l)-z"2.[(l-Aa)  +  (l+Aa)ic]|  __  ,=  0. 


z=-l 

<=1 


Thus,  approximation  (**.7)  is  unstable  due  to  the  violation  of  the  necessary 
stability  condition  (3-U)  at  the  point  z  =  -1. 

Recalling  the  summary  result  which  follows  Example  h.6,  we  find  that 
approximation  (U.7)  satisfies  all  the  required  hypothesis  except  for  condition 
(U.6d).  Indeed,  we  have 


(b.7c 


P(z=-l,tc=l)  =  1-z  ^«Aa(ic-ic  1)-z2*[l-j^(<-l)2(l-ic  ^)2] 


z=-l 

k=1 


=  0. 


We  remark  that  according  to  Lemma  3.2,  z  =  -1  is  the  only  possible  point  for  the 
thrae-lcvel  basic  scheme  (l*.7a)  to  violate  Assumption  V,  as  we  have  indeed  found 
in  (b.7c) • 

We  close  this  chapter  by  considering  difference  approximations  to  the  two 


* 


apace  dimensional  problem 


I 
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3u(x ,y ,t ) / 3t  =  a3u(x,y,t)/3x  +  b3u(x  ,y  ,t )/3y  ,  a>  0  ;  u(x,y,0)  =  f(x,y), 

in  the  quarter  space  x  >,  0 ,  t  z  0 ,  -«  <  y  <  ® .  The  analysis  of  such  initial¬ 
boundary  problems  in  both  the  differential  case  (see  [9])  and  the  difference 
case  (see  for  example  [l])  can  be  carried  out  by  Fourier  transforming  vitn 
respect  to  the  variable  y;  thus  obtaining  a  one  space  dimensional  problem  of 
the  type  analyzed  in  the  previous  chapters.  To  be  more  precise,  let  Ax,  Ay 

be  the  spatial  mesh  width  such  that  x  =  At/Ax  =  constant, 

x 

X  =  At/Ay  =  constant  and  denote  by  v  (t)  =  v(vAx,£Ay,t)  the  approximated 

y  ^  c 

grid  function.  Then,  Fourier  transforming  in  the  y-direction  (with  dual 
variable  n)  and  Fourier-Laplace  transforming  with  respect  to  the  tine-variable 
t  leads  one  to  search  for  normal  modes  of  the  form  v  (t)  -  znpc^ei^n.  Upon 
substituting  such  modes  as  a  grid  solution  for  a  given  difference  scheme, 
we  obtain  the  corresponding  characteristic  equation.  If  in  particular  the  scheme 
is  the  one  generates  the  translatory  boundary  conditions  considered,  we  obtain 
the  associated  boundary- function  which  determines  the  stability  properties  of 
these  conditions.  Both  the  characteristic  and  boundary  functions  involved  in 
the  two  space  dimensional  case  are  dependent  on  the  extra  parameter  n,  and 
our  former  results  are  still  valid  in  this  case  since  all  estimates  made  are 
unifonu..in  n  (see  [6], [8]).  The  only  exception  is  that  of  Lemma  3.1,  according 
to  which  the  possibility  of  k  =  1  to  serve  as  an  inner  solution  corresponding 
to  z  =  1,  is  excluded  in  the  outflow  case  (a>0).  The  proof  of  Lemma  3.1  is 
based  on  the  consistency  condition,  so  it3  validity  in  the  two  space  dimensional 
case  is  restricted  to  the  single  point  n  *  0.  Therefore,  for  the  result  of 


I 
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Lemma  3- 1  to  be  valid  independently  of  the  extra  parameter  n,  the  additional 
requirement 

(**.8)  RQ(z=i,ic=i,n)  ¥  o  ,  o  <  |n|  $ 


must  be  fulfilled.  Thus,  to  apply  our  scheme-independent  stability  criteria  for 
a  two  space  dimensional  approximation,  we  first  determine  the  boundary  stability 
properties  by  employing  the  associated  boundary-function.  Then  it  remains  to 
check  whether  the  approximation  meets  the  additional  assumptions  V  and  VI,  and 
whether  condition  (4.8)  is  fulfilled. 

We  note  that  when  verifying  the  validity  of  Assumption  V  in  the  two  space 
dimensional  case,  one  may  no  longer  use  Lemma  3.2  in  which  conditions  for 
the  validity  of  Assumption  V  for  tvc- and  three-leveled  schemes,  are  discussed. 
Indeed,  the  lemma  follows  from  Lemma  2.2  whose  proof  is  based  on  the  zero  order 
accuracy  condition,  i.e.,  when  dealing  with  the  boundary  conditions  we  have 
R0(z=l,ie=»l,n=0)  ~  0.  Thus,  the  conclusions  of  Lemma  2.2  and  Lemma  3.2  hold 
only  in  the  neighbourhood  of  n  *  0  and  are  not  necessarily  valid  for  all  n, 

0  $  1  n  1  $  w. 

EXAMPLE  4.7.  Let  the  boundary  conditions  be  generated  by  the  right-sided  explicit 
Euler  scheme, i.e., 


v  (t+At) 

W  9  S 


yCv„.c*i<t)-v<(tn 


i 


I 


( 4.9a) 


Vi  =  0,-1,. . .  ,-r+l,  -«  <  ;  <  00 . 


I 
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"he  boundary- function  associated  with  (U.9a)  is  given  by 


(4.9b)  R0(z,tc,n)  =  1  -z  *[l+X  a*(<-l)+X  b-(e  -l)], 

u  x  y 


r.:1  by  oquatin**  to  vo  r^t  that  its  z  -solutions  satisfy 


|  z(  n )  I  2=  ( l->  a-X  t>+X  a*cos£+X  b-cosn)2  *  (  X  a-sir.^+X  bv.im)"  < 

1  x  y  x  y  x  y 

(4.9c) 

<  (|l->xa-^ybl  +  Uxa+Xybj)2  ,  0  <  |c|  f  it,  Of  |  n  |  f  it. 


Then,  upon  imposing  the  CFL  condition 


(U.9d) 


0  <  X  a  +  X  b  f  1 
x  y 


we  see  that  the  explicit  and  hence  (by  Remark  3.3)  solvable  boundary  conditions 
(4.9a)  are  of  dissipative  type.  By  (4.9d)  it  follows  that  in  the  outflow  case 
(a>0)  we  have  X^b  f  1,  hence 


|  zU=l,n)  |  =  (l-X  b+X  b*cosn)  +  (X  b-sinn;  < 1 1— X  bj+|X  b|=*  1, 

y  y  y  y  y 


( 4 . 9e ) 


0  <  [h  |  ?  tt  , 


so  condition  (4.8)  is  fulfilled.  Consistency  implies  that  z(ic»l,n*0)  »  1  (see 
(4.9e))  and  together  with  (4.9e)  we  finally  get 


I 
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A(z=ei<P,n)  i  |  Rq(  z=e'l'v,<=l  ,n )  j  >  0  ,  0  <  |cp|  *  it  ,  0  (  |  n  (  £  * 


i(p 


i.e.,  Assumption  V  holds.  Therefore  we  may  apply  Theorem  3-4  concluding  that 
if  the  CFL  condition  (4.9d)  is  fulfilled then  the  outflow  boundary  conditions 
(4.9a)  always  maintain  stability  independently  of  the  interior  scheme. 

Boundary  conditions  ( U . 9a. )  are  generated  by  an  obvious  extension  of  the  usual 
one  space  dimensional  right-sided  Kuler  scheme  discussed  in  Exanrale  4.4.  Another 
possible  extension  is  given  by  (see  example  (2.5)  in  [l]) 


(4.10a) 


u  *  0,-1, . . , ,-r+l,  -®  <  C  < 


These  boundary  conditions  are  unstable  in  the  sense  that  the  z-solutions  of  the 
associated  boundary- function 


(4.10b) 


RQ(z,x,n) 


1-z  *[X  a(ic-l)  +  iX  b.sinnl, 

x  y 


satisfy 


(4.10c)  |z(x=e  ^,n)|^  3  (l-X^a+X^a^cosO^  +  (X  a*sin£  +  X  b*sinn)^> 

x  y 


x  x 


hence  jz^se1  , n )  |  1 2  g  4  »  (Xyb*alnn)2  >  1,  O'  <  |n|  S  *  .  Thus,  boundary 


(l)  A  further  restriction  on  the  time  step  At  may  arise  from  the  stability 
requirement  made  in  Assumption  III  (see  Remark  4.1), 
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conditions  (l+.lOa)  fail  to  satisfy  the  von  Neumann  condition  for  ail  n,  Of|n|jit 
and  our  scheme-independent  stability  criteria  are  inapplicable  in  this  case. 

The  question  of  stability  is,  in  this  case,  dependent  on  the  basic  scheme  utilized. 

A  further  possibility  to  extend  the  one  space  dimensional  Euler  scheme  (4.4a) 
is  considered  in  the  following  example. 

EXAMPLE  4.8.  (example  (2.6)  in  Cl]).  Let  the  boundary  conditions  be  of  the  form 


v  (t+At)  =  v  (t)  +  0.5‘X  a*  C(v  .  ,(t)+v  ,  (t))-(v  .(t)+v  (t)]  + 

u,G  x  y+l,C+l  u+l,C  P.C+l  u,C 


(4.11a) 


*  °'5' V  ' 


JJ  —  0,-1  ,  •  .  .  1*+1  ,  — 00  <  £  <*  00^. 


To  simplify  the  computations  we  shall  consider  the  case  whe'ffe  A^a  *  A^b.  ’fhe 
boundary-function  associated  with  (4.11a)  is  given  by 


(4.11b) 


R0(z,x,n) 


l-8“1*Xxa‘(«fein-l) 


and  by  equating  to  zero  we  get  that  its  z-solutions  satisfy 


(4.11c)  | z( K=e*^ ,n) | 2  =  [l-A  a+A  a*cos(^+h)]^  +  [A  a*sin(£*n)]2  J 


x  x 


f  ( 1 1— A^a.  |  +|Axa|)2  ,  0  €  |V|,|n|  *  *  • 


Then,  upon  imposing  the  CFL  condition 


I 
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f'-.lld)  0<Xa=Xbsl 

x  y 

we  see  that  the  explicit  and  hence  (by  Rorark  3-3)  solvable  boundary  conditions 

{ h.  1  la)  rrt.i r  fy  the  vnn  ’•l- unarm  condition.  ?-y  (I*. 11c)  end  (h.lld)  we  have 


'  ■  -1 , n ) | ^  =  ( 1 -X  a+X  a-cosn)^  +  (X  a*sinn)^  <  ( 1 1— X  a|  +  |x  a| )2  =  1, 
xx  x  xx 

(-. Lie) 

o  <  | n |  s  7T, 

and  hence  condition  (4.8)  is  fulfilled.  Consistency  implies  that 
z(tc=l,n=0)  =  1  (see  (4. lie))  and  together  with  (4. lie)  we  finally  get 

A(z=ell45,n)  i  |R0(z=eltp,K=l,n)|  >  0  ,  0  <  |<p|  s  w ,  0  f  | n |  (  if , 


i.e..  Assumption  v  holds.  Therefore  we  may  apply  Theorem  3.3  concluding  that 
if  the  CFL  condition  (4. lid)  is  fulfilled[ ^  then  the  outflow  boundary  conditions 
(4.11a)  in  conjunction  with  any  dissipative  basic  scheme  constitute  together  a 
stable  approximation. 

In  the  manner  of  the  last  two  examples,  one  may  consider  various  two  space 
dimensional  boundary  treatments  which  extend  the  corresponding  one  space 
dimensional  ones.  We  choose  to  consider  an  example  which  is  based  on  extrapolation. 
EXAMPLE  4.9.  (example  (2.4)  in  [l]).  Let  the  boundary  conditions  be  determined 

(l)  A  further  restriction  on  the  time  step  At  may  arise  from  the  stability 
requirement  made  in  Assumption  III  (see  Remark  4.1). 
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by  oblique  extrapolation  along  the  characteristic  plane  (see  Example  U.l),  i.e.. 


m 

(  4. 12a'  v  ( t+At )  =  :‘T)(-1)J+1ViiA,  Ct-(j-l)At],  p=0,-l,. ..  ,-r+l, -»<;«?>. 


W+J.C+J 


Phe  boundary- function  associated  with  (U.12a)  is  given  by 


( 12b) 


^0(z,ie,n)  =  1 


-  II  (")(-DJ*VJ 

J-l  J 


-JJ  ij«  -  /,  -1  in 


-  /,  -J.  xr|\ 

e  =  (1-z  <e  )  , 


and  by  equating  to  zero  we  get  that  its  z-solutions  satisfy 


{U.  l:V; 


zU=eU,n)|  *  |eUein|  »  1,  Os  U|  ,  |  n  |  s  it 


Thus,  the  explicit  and  hence  (by  Remark  3.3)  solvable  boundary  conditions  (U.12a) 
are  of  non-dissipative  type  ao  they  satis*y  the  von  Neumann  condition.  In 
addition  we  have 

in  m 

R0(z«l,K»l,n)  ■  (l-e  )  i  0  ,  0  <  j n |  *  ir, 

i.e.,  condition  (U.8)  is  fulfilled.  Thus  to  assure  stability  via  Theorem  3.3S  it 
then  remains  to  verify  Assumption  V.  Since  by  (U.12b)  we  have 


R0(z=eio,|c=l,n)|^n  -  (l-e“1<peln)V  n  -  0, 


the  assumption  is  reduced  in  this  case  to  the  requirement 
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(U.l2d)  A(  z=e1(P,n )  =  |P(z=eliP,ic=l,n)  |  >  0,  0  <  |<p|  s  w,  0  <  |n|  £ 

P(z,<,n)  denoting  as  usual  the  corresponding  characteristic  function  associated 
with  the  basic  scheme. 

Applying  Theorem  3.3  we  conclude  that  the  extrapolatory  outflow  boundary 
conditions  (4.12a) ,  when  complementing  any  dissipative  basic  scheme  whose 
characteristic  function  P(z,K,n)  satisfies  P( z“e1(^,ie*l »n)  4  0,  0<|<p|<*, 

O'  I  n  |V*  ,  constitute  a  stable  approximation. 


\ 

\ 


-  g 


3  - 


REFERENCES 


Cl]  Abarbanei  and  D.  Gottlieb,  Stability  of  two  dimensional  Initial  boundary 

value  problems  using  leap  Frog  type  schemes,  Dept,  of  Math.  Sciences, 

Tel -Aviv  'Jniv.  P.ep.  ,  1977. 

1.9]  T.  Flvins  and  A.  Sundstrom,  Computationally  efficient  schemes  and  boundary- 
conditions  for  a  fine-mesh  barotropic  model  based  on  the  shallow-water 
equations,  Tellus,  v.  25,  1973,  pp.  132-156. 

[R]  M.  Goldberg,  On  a  boundary  extrapolation  theorem  by  Kreiss,  Math.  Comp. , 
v.  ^1,  1977,  pp.'  U^9-UT7. 

[(‘]  M.  Goldberg  and  E.  Tadmor,  Scheme-independent  stability  criteria  for 

difference  approximations  of  hyperbolic  initial-boundary  value  problems,  I, 

Math,  comp.,  v.  22,  1978,  pp.  1097-1107. 

[5]  3.  Gustafsson,  On  the  convergence  rate  for  difference  approximations  to  mixed 
Initial  boundary  value  problems ,  Math.  Comp.;  v.  29.  1975.  op.  996-lQo] 

[6]  B.  Gustafsson,  H.O.  Kreiss  and  A.  Sundstrttm,  Stability  theory  of  difference 
approximations  for  mixed  initial-boundary  value  problems.  II,  Math.  Comp. , 
v.  _2&,  l972,  pp  6U9-686. 

[7.1  H.O.  Kreiss,  Difference  approximations  for  hyperbolic  differential  equations. 

Proc.  Sympos.  Univ.  of  Maryland,  Maryland,  Acad.  Press,  New  York,  1966,  PP- 51-53. 

[3]  H.O.  Kreiss,  Stability  theory  for  difference  approximations  of  mixed  initial 
boundary  value  problem.  I.  Math.  Comp.,  v.  22,  1968,  pp.  703-71*+,  MR  39^2355. 

[9]  H.O.  Kreiss,  Initial  boundary  value  problems  for  hyperbolic  systems.  Comm. 

Pure.  Appl.  Math.,  v.  23.  1970.  pp.  277-298. 

ClO]  H.O.  Kreiss  and  E.  Lundqyist,  On  difference  approximations  with  wrong  boundary 
values .  Math.  Comp.,  v.  _22_,  1968,  pp.  1-12. 

Cll]  H.O.  Kreiss  and  J.  Oliger,  Methods  for  the  approximate  solution  of  time 
dependent  problems.  GARP  publication  series  No.  10,  1973. 

[12]  S.  Osher,  Stability  of  parabolic  difference  approximations  to  certain  mixed 
initial  boundary  value  problems.  Math.  Comp. .  v.  26,  1972.  pp.  13-39. 

[13]  R.D.  Richtmyer  and  K.W.  Morton,  Difference  Methods  for  Initial  Value  Problems, 

Interscience,  New  York,  1967-  " 


[l*+]  G.  SkBllermo,  How  the  boundary  conditions  affect  the  stability  and  accuracy  of 


-  89  - 


some  implicit  methods  for  hyperbolic  equations,  Rep.  So.  62,  1975,  Dept,  c 
Comp.  Sci.  Uppsala  Univ. ,  Uppsala,  Sweden. 

[l5l  9.  SkttllernD,  Error  analysis  for  the  mixed  initial  boundary  value  problem 
for  hyperbclic~~equations ,  Rep.  No.  ^3.  1975.  Dept,  of  Comp.  Sci..  Uppsala 
'Jniv.  ,  Uppsala,  Sweden. 

[ 16]  9.  Strang,  Implicit  difference  methods  for  initial -boundary  value  problems, 
J.  Math.  Aral.  Appl. .  v.  l6.  1966.  pp.  188-198. 


END 


DATE 

FILMED 


DTIC 


